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ABSTRACT 

In  response  to  the  transition  by  the  United  States  Nuclear  Regulatory  Commission  (NRC)  to  a 
risk-informed,  performance-based  fire  protection  rulemaking  standard,  Fire  Probabilistic  Risk  Assessment 
(PRA)  methods  have  been  improved,  particularly  in  the  areas  of  advanced  fire  modeling  and  computational 
methods.  As  the  methods  for  the  quantification  of  fire  risk  are  improved,  the  methods  for  the  quantification 
of  the  uncertainties  must  also  be  improved.  In  order  to  gain  a  more  meaningful  insight  into  the  methods 
currently  in  practice,  it  was  decided  that  a  scenario  incorporating  the  various  elements  of  uncertainty 
specific  to  a  fire  PRA  would  be  analyzed. 

The  NRC  has  validated  and  verified  five  fire  models  to  simulate  the  effects  of  fire  growth  and 
propagation  in  nuclear  power  plants.  Although  these  models  cover  a  wide  range  of  sophistication, 
epistemic  uncertainties  resulting  from  the  assumptions  and  approximations  used  within  the  model  are 
always  present.  The  uncertainty  of  a  model  prediction  is  not  only  dependent  on  the  uncertainties  of  the 
model  itself,  but  also  on  how  the  uncertainties  in  input  parameters  are  propagated  throughout  the  model. 
Inputs  to  deterministic  lire  models  are  often  not  precise  values,  but  instead  follow  statistical  distributions. 

The  fundamental  motivation  for  assessing  model  and  parameter  uncertainties  is  to  combine  the 
results  in  an  effort  to  calculate  a  cumulative  probability  of  exceeding  a  given  threshold.  This  threshold  can 
be  for  equipment  damage,  time  to  alarm,  habitability  of  spaces,  etc. 

Fire  growth  and  propagation  is  not  the  only  source  of  uncertainty  present  in  a  fire-induced 
accident  scenario.  Statistical  models  are  necessary  to  develop  estimates  of  lire  ignition  frequency  and  the 
probability  that  a  lire  will  be  suppressed.  Human  Reliability  Analysis  (HRA)  is  performed  to  determine  the 
probability  that  operators  will  correctly  perform  manual  actions  even  with  the  additional  complications  of  a 
fire  present. 

Fire  induced  Main  Control  Room  (MCR)  abandonment  scenarios  are  a  significant  contributor  to 
the  total  Core  Damage  Frequency  (CDF)  estimate  of  many  operating  nuclear  power  plants.  Many  of  the 
resources  spent  on  fire  PRA  are  devoted  to  quantification  of  the  probability  that  a  fire  will  force  operators 
to  abandon  the  MCR  and  take  actions  from  a  remote  location.  However,  many  current  PRA  practitioners 
feel  that  effect  of  MCR  fires  have  been  overstated. 

This  report  details  the  simultaneous  application  of  state-of-the-art  model  and  parameter 
uncertainty  techniques  to  develop  a  defensible  distribution  of  the  probability  of  a  forced  MCR 
abandonment  caused  by  a  fire  within  a  MCR  benchboard.  These  results  are  combined  with  the  other 
elements  of  uncertainty  present  in  a  fire-induced  MCR  abandonment  scenario  to  develop  a  CDF 
distribution  that  takes  into  account  the  interdependencies  between  the  factors.  In  addition,  the  input  factors 
having  the  strongest  influence  on  the  final  results  are  identified  so  that  operators,  regulators,  and 
researchers  can  focus  their  efforts  to  mitigate  the  effects  of  this  class  of  fire-induced  accident  scenario. 

Thesis  Supervisor:  George  E.  Apostolakis 

Professor  of  Nuclear  Science  and  Engineering  and  Engineering  Systems 


DISCLAIMER 


The  views  expressed  in  this  doeument  are  those  of  the  author  and  do  not 
neeessarily  represent  the  position  of  the  Department  of  Defense  or  any  of  its  components, 
including  the  Department  of  the  Navy. 


ACKNOWLEDGEMENTS 


I  would  like  to  thank  my  advisor,  Professor  George  Apostolakis,  for  his  guidanee 
and  patienee  throughout  the  duration  of  this  project. 

I  would  also  like  to  thank  all  of  the  NRC  staff  who  have  offered  valuable 
technical  advice  during  the  performance  of  this  work.  Specifically,  I  would  like  to  thank 
Don  Helton  and  Nathan  Siu  for  their  helpful  suggestions  and  insights. 

Furthermore,  I  am  especially  grateful  to  my  fellow  MIT  student  Dustin 
Langewisch  for  his  invaluable  advice  and  assistance  with  some  of  the  more 
computationally  intensive  portions  of  this  study. 


TABLE  OF  CONTENTS 


I.  EXECUTIVE  SUMMARY . 1 

II.  SCENARIO  IDENTIEICATION . 6 

III.  EIRE  IGNITION  EREQUENCY . 16 

IV.  EIRE  SCENARIO  MODELING . 21 

V.  MODEL  UNCERTAINTY . 30 

VI.  PARAMETER  UNCERTAINTY . 42 

VIE  COMBINED  MODEL  AND  PARAMETER  UNCERTAINTY . 51 

VIII.  EIRE  SUPPRESSION  ANALYSIS . 55 

IX.  HUMAN  RELIABILITY  ANALYSIS . 60 

X.  SENSITIVITY  ANALYSIS  AND  CONCLUSTIONS . 64 

REEERENCES . 73 

APPENDIX  A:  EMSNL  INPUT  PILE . 76 

APPENDIX  B;  EMSNL  21  INPUT  TO  CPAST-SCREEN  VIEW . 77 

APPENDIX  C:  EMSNL  THERMAL  PROPERTIES  INPUT  PILE . 82 

APPENDIX  D:  WINBUGS  INPUT  PILE  POR  HGL  DATA  (BE3  ONLY) . 83 

APPENDIX  E:  WINBUGS  INPUT  PILE  POR  HGL  DATA  (BE3  +  EMSNL 

21/22) . 84 

APPENDIX  P:  EMSNL  21  INPUT  TO  PPS . 85 

APPENDIX  G:  STUDENT  BIOGRAPHY . 88 


IV 


LIST  OF  FIGURES 


Figure  2-1:  RCP  Seal  LOCA  Sequence  of  Events . 7 

Figure  2-2:  CPSWPR  Sequence  of  Events . 10 

Figure  2-3:  Event  Tree  T3-  Turbine  Trip  with  MFW  Available .  1 1 

Figure  2-4:  Event  Tree  S2-Small  LOCA .  12 

Figure  2-5:  Fire  Effects  on  Internal  Sequence  of  Events . 12 

Figure  3-1:  MCR  Fire  Ignition  Frequency  Given  inNUREG-1150 . 17 

Figure  3-2:  MCR  Fire  Ignition  Frequency  Given  in  EPRI  1016735  .  18 

Figure  3-3:  Area  Ratio  of  Benchboard  1-1  to  total  MCR  Cabinet  Area . 19 

Figure  3-4:  Frequency  of  Fires  in  Benchboard  1-1 . 20 

Figure  4-1:  FMSNL  21  Heat  Release  Rate . 25 

Figure  4-2:  CFAST  HGL  Height . 26 

Figure  4-3:  CFAST  Optical  Density . 27 

Figure  4-4:  CFAST  Heat  Flux  to  Operator . 28 

Figure  4-5:  CFAST  HGL  Temperature . 28 

Figure  5-1:  Comparison  of  Model  Prediction  to  Experimentally  Determined  HGL  Temperature . 31 

Figure  5-2:  CFAST  HGL  Temperature  Prediction  vs.  BE3  Experimental  Data . 35 

Table  5-1:  UMD  BE3  WinBUGS  Data . 35 

Figure  5-4:  CFAST  HGL  Temperature  Prediction  vs.  BE3  and  FMSNL  21/22  Experimental  Data . 37 

Figure  6-1 :  Comparison  of  CFAST  Output  to  MQH  Prediction . 43 

Figure  6-2:  Comparison  of  HRR  Distributions  Given  Variations  in  Growth  Time . 46 

Figure  6-3:  NUREG-6850  Appendix  G  HRR  for  Cabinets  with  Qualified  Cable . 47 

Table  6-4:  Random  Variables  Used . 47 

Figure  6-4:  Convergence  of  Probabilistic  Fire  Simulator  Results . 48 

Figure  6-5:  Peak  Predicted  HGL  Temperature  Distribution  Prediction . 49 

Figure  6-6:  Cumulative  Distribution  Function  of  Peak  Predicted  HGL  Temperature . 50 

Figure  7-1:  Comparison  of  Model  Uncertainty  Techniques  by  Mean  Peak  HGL  Temperature . 52 

Figure  7-2:  CDF  of  Forced  MCR  Abandonment . 53 


V 


Figure  7-  3:  Comparison  of  Model  Uncertainty  Techniques  by  Probability  of  Exceedance . 54 

Figure  8-1:  Fire  Suppression  Event  Tree . 56 

Figure  8-2:  FMSNL  Time  Available  for  Suppression . 57 

Figure  8-3:  Probability  of  Non-Suppression . 57 

Figure  8-4:  MCR  Abandonment  Probability  Given  in  NUREG-1 150 . 59 

Figure  9-1 :  Scoping  FIRA  Analysis  for  MCR  Abandonment  Scenario . 61 

Figure  9-2:  Probability  of  Successful  Operator  Action  from  the  RSP  Given  in  NUREG-1 150 . 62 

Figure  10-1:  Comparison  of  CDF  Distributions . 65 

Table  10-1:  Comparison  of  CDF  by  Method . 65 

Table  10-2:  Comparison  of  MCR  Abandonment  Probability  by  Method . 66 

Figure  10-2:  Comparison  of  MCR  Abandonment  Distributions . 67 

Figure  10-3:  Rank  Order  Correlation  Coefficient  for  Peak  FIGL  Temperature . 68 

Figure  10-4:  FMSNL  21  FIGL  Temperatures  as  a  Function  of  Ventilation  Rate . 69 

Figure  10-5:  Effect  of  Doubling  Fire  Suppression  Rate . 70 


VI 


LIST  OF  ACRONYMS 


ATHEANA 

AUX  BLDG 

BE3 

CCW 

CDF 

CFAST 

CPSWPR 

CV/T 

EPRI 

ESWGR 

FDS 

FDTs 

FIVE-Rev1 

FMAG 

FMSNL 

HEP 

HFE 

HGL 

HPI 

HRA 

HRR 

LOCA 

MGR 

MIT 

NIST 

NRG 

PDF 

PFS 

PORV 

PRA 

RCC 

RCP 

RES 

RSP 

SA 

UA 

UMD 

WINBUGS 


A  TECHNIQUE  FOR  HUMAN  EVENT  ANALYSIS 
AUXILIARY  BUILDING 
BENCHMARKING  EXERCISE  #3 
COMPONENT  COOLING  WATER 
CORE  DAMAGE  FREQUENCY 

CONSOLIDATED  MODEL  OF  FIRE  SMOKE  AND  TRANSPORT 

CHARGING  PUMP  SERVICE  WATER  PUMP  ROOM 

CABLE  VAULT/TUNNEL 

ELECTRIC  POWER  RESEARCH  INSTITUE 

EMERGENCY  SWITCHGEAR  ROOM 

FIRE  DYNAMICS  SIMULATOR 

FIRE  DYNAMICS  TOOLS 

FIRE-INDUCED  VULNERABILITY  EVALUATION,  REVISION  1 

NUCLEAR  POWER  PLANT  FIRE  MODELING  APPLICATION  GUIDE 

FACTORY  MUTUAL/SANDIA  NATIONAL  LABORATORY 

HUMAN  ERROR  PROBABILITY 

HUMAN  FAILURE  EVENT 

HOT  GAS  LAYER 

HIGH  PRESSURE  INJECTION 

HUMAN  RELIABILITY  ANALYSIS 

HEAT  RELEASE  RATE 

LOSS  OF  COOLANT  ACCIDENT 

MAIN  CONTROL  ROOM 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
NATIONAL  INSTITUE  OF  STANDARDS  AND  TECHNOLOGY 
UNITED  STATES  NUCLEAR  REGULATORY  COMMISSION 
PROBABILITY  DENSITY  FUNCTION 
PROBABILISTIC  FIRE  SIMULATOR 
POWER  OPERATED  RELIEF  VALVE 
PROBABILISTIC  RISK  ASSESSMENT 
RANK  ORDER  CORRELATION  COEFFICIENT 
REACTOR  COOLANT  PUMP 

U.S.  NRC  OFFICE  OF  NUCLEAR  REGULATORY  RESEARCH 

REMOTE  SHUTDOWN  PANEL 

SENSITIVITY  ANALYSIS 

UNCERTAINTY  ANALYSIS 

UNIVERSITY  OF  MARYLAND 

WINDOWS  BAYESIAN  INFERENCE  USING  GIBBS  SAMPLING 


I,  EXECUTIVE  SUMMARY 


In  response  to  the  transition  by  the  United  States  Nuelear  Regulatory  Commission 
(NRC)  to  a  risk-informed,  performanee-based  fire  protection  rulemaking  standard  [1], 

Fire  Probabilistic  Risk  Assessment  (PRA)  methods  have  been  improved  [2],  particularly 
in  the  areas  of  advanced  fire  modeling  and  computational  methods  [3].  As  the  methods 
for  the  quantification  of  fire  risk  are  improved,  the  methods  for  the  quantification  of  the 
uncertainties  must  also  be  improved.  In  order  to  gain  a  more  meaningful  insight  into  the 
methods  currently  in  practice,  it  was  decided  that  a  scenario  incorporating  the  various 
elements  of  uncertainty  specific  to  a  fire  PRA  would  be  analyzed. 

Fire-induced  Main  Control  Room  (MCR)  abandonment  scenarios  are  a  significant 
contributor  to  the  total  Core  Damage  Frequency  (CDF)  estimate  of  many  operating 
nuclear  power  plants  [4].  Many  of  the  resources  spent  on  fire  PRA  are  devoted  to 
quantifying  the  probability  that  a  fire  will  force  operators  to  abandon  the  MCR  and  take 
actions  from  a  remote  location.  However,  many  current  PRA  practitioners  [3]  feel  that 
the  effects  of  MCR  fires  have  been  overstated.  This  thesis  demonstrates  the  application  of 
state-of-the-art  techniques  for  analyzing  the  uncertainty  and  sensitivity  of  a  fire-induced 
MCR  abandonment  scenario. 

The  NRC  has  validated  and  verified  five  fire  models  to  simulate  the  effects  of  fire 
growth  and  propagation  in  nuclear  power  plants  [5].  Although  these  models  cover  a  wide 
range  of  sophistication,  epistemic  uncertainties  resulting  from  the  assumptions  and 
approximations  used  within  the  model  are  always  present. 

For  our  scenario,  the  Consolidated  Model  of  Fire  Growth  and  Smoke  Transport 
(CFAST)  [6]  is  used  to  predict  the  evolution  of  environmental  conditions  after  the 


1 


ignition  of  a  MCR  fire  (Section  IV).  CFAST  was  chosen  because  adequate  and 
computationally  inexpensive  results  can  be  obtained  for  simple  configurations  like  ours 
[5].  The  primary  simplification  inherent  to  CFAST  is  the  assumption  that  each 
compartment  can  be  subdivided  into  two  zones  that  are  uniform  in  temperature  and 
species  concentration.  Choosing  a  so-called  zone  model,  such  as  CFAST,  allows  for 
much  larger  Monte  Carlo  samples  to  be  reasonably  achieved  in  determining  model  input 
parameter  uncertainties  (Section  VI)  and  sensitivity  analyses  (Section  X). 

The  upper  zone  in  a  zone  model  is  referred  to  as  the  Hot  Gas  Layer  (HGL).  Of 
the  MCR  abandonment  criteria  [2],  the  results  of  this  study  indicate  that  the  peak  HGL 
temperature  reached  is  the  limiting  factor  in  predicting  forced  MCR  abandonment.  For 
our  scenario,  the  evolution  of  the  environmental  conditions  predicted  by  CFAST  reveal 
that  the  HGL  layer  height  will  descend  rapidly  after  fire  ignition,  while  the  HGL 
temperature  will  take  several  additional  minutes  to  reach  a  value  that  would  force 
evacuation. 

The  general  method  of  evaluating  model  uncertainty  is  through  the  comparison  of 
model  data  with  that  of  actual  experiments.  A  Bayesian  Framework  for  Model 
Uncertainty  Considerations  in  Fire  Simulation  Codes  [7],  proposed  by  the  University  of 
Maryland  (UMD),  was  first  conducted  by  comparing  CFAST  output  to  experimental 
results  contained  in  NUREG-1824’s  Benchmarking  Exercise  Three  (BE3)  [5].  Through 
cooperation  with  the  United  States  Nuclear  Regulatory  Commission  (NRC)  and  the 
Massachusetts  Institute  of  Technology  (MIT),  the  UMD  method  was  extended  to  include 
data  from  the  Factory  Mutual/Sandia  National  Laboratory  (FMSNL)  21  and  22  tests 
conducted  as  part  of  NUREG/CR-4527,  An  Experimental  Investigation  of  Internally 
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Ignited  Fires  in  Nuclear  Power  Plant  Control  Cabinets  Part  II:  Room  Effects  Tests  [8]. 
The  results  from  the  UMD  method  are  then  compared  to  the  method  presented  in 
NUREG-1934,  Nuclear  Power  Plant  Fire  Modeling  Application  Guide  (FMAG)  [3]. 

Comparison  of  test  data  from  [8,  9]  to  model  predictions  shows  that  CFAST 
consistently  over-predicts  HGL  temperature.  Therefore,  the  most  conservative  method  of 
analyzing  our  scenario  would  be  to  neglect  model  uncertainty  altogether  and  use  the 
values  predicted  by  CFAST.  Less  conservative  results  can  be  obtained  by  using  the 
method  presented  in  the  FMAG,  with  the  UMD  method  yielding  the  least  conservative 
results. 

The  uncertainty  of  a  model  prediction  is  not  only  dependent  on  the  uncertainties 
of  the  model  itself,  but  also  on  how  the  uncertainties  in  input  parameters  are  propagated 
throughout  the  model.  Inputs  to  deterministic  fire  models  are  often  not  precise  values,  but 
instead  follow  statistical  distributions.  Due  to  the  complexity  and  non-linear  nature  of 
our  fire  model,  empirical  methods  to  estimate  uncertainty  propagation  do  not  yield 
sufficiently  refined  results  when  multiple  input  parameters  are  allowed  to  vary.  In  order 
to  more  adequately  assess  the  distributions  of  fire  model  output  variables,  Monte  Carlo 
simulations  have  been  coupled  to  our  fire  model,  CFAST,  through  a  tool  called 
Probabilistic  Fire  Simulator  (PFS)  [10]. 

The  fundamental  motivation  for  assessing  model  and  parameter  uncertainties  is  to 
combine  the  results  in  an  effort  to  calculate  a  cumulative  probability  of  exceeding  a  given 
threshold.  This  threshold  can  be  for  equipment  damage,  time  to  alarm,  or,  for  our 
scenario,  the  habitability  of  a  MCR  due  to  HGL  temperature. 
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Combining  current  model  uneertainty  methods  (FMAG,  UMD)  with  parameter 
uncertainties  (PFS)  results  in  a  reduetion  by  a  faetor  of  approximately  20  in  the  mean 
value  of  the  probability  of  foreed  MCR  abandonment  ealeulated  in  NUREG-1 150. 

In  evaluating  the  eombined  model  and  parameter  uneertainties  present  in  our 
scenario,  the  goal  was  to  develop  an  expression  for  the  probability  that  a  fire  in  a 
benehboard  would  foree  operator  abandonment  of  the  MCR  if  no  suppression  efforts 
were  made.  However,  fire  growth  and  propagation  is  not  the  only  souree  of  uneertainty 
present  in  a  fire-indueed  aceident  scenario.  Statistical  models  are  neeessary  to  develop 
estimates  of  fire  ignition  frequeney  [11]  (Seetion  III)  and  the  probability  that  a  fire  will 
be  suppressed  [2]  (Seetion  VIII).  Human  Reliability  Analysis  (HRA)  [12]  (Seetion  IX)  is 
performed  to  determine  the  probability  that  operators  will  eorreetly  perform  manual 
actions  even  with  the  additional  eomplications  stemming  from  the  presenee  of  a  fire. 

The  elements  of  uneertainty  present  in  a  fire-induced  MCR  abandonment  scenario 
are  combined  to  develop  a  CDF  distribution  that  takes  into  aeeount  the  interdependeneies 
between  the  faetors.  The  eurrent  methods  used  in  this  study  show  that  a  reduetion  by  a 
faetor  of  approximately  two  in  the  mean  value  of  the  total  CDF  calculated  in  NUREG- 
1 150  is  expeeted.  Although  this  value  is  lower  than  previously  assessed  [4],  it  is  still  a 
signifieant  eontributor  to  total  CDF  and  would  not  be  eliminated  in  the  sereening  proeess 
of  a  fire  PRA. 

An  Uneertainty  Analysis  (UA)  is  ineomplete  without  a  diseussion  of  whieh  input 
faetors  have  the  strongest  infiuenee  on  the  results.  Sensitivity  Analysis  (SA)  is  defined 
by  Saltelli,  et  al.  [13]  as  “[t]he  study  of  how  uncertainty  in  the  output  of  a  model 
(numerical  or  otherwise)  can  be  apportioned  to  different  sources  of  uncertainty  in  the 
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model  input.  ”  There  are  many  SA  methods  available  [14],  and  choosing  the  one  that  is 
most  appropriate  is  important  to  yield  meaningful  results.  In  order  to  determine  which 
input  parameters  have  the  strongest  influence  on  the  results  of  our  fire  model,  PFS  uses 
the  Spearman  Rank-order  Correlation  Coefficient  (RCC)  [15]  to  assess  the  sensitivity  of 
an  output  value  to  an  input  parameter.  RCC  is  independent  of  the  distribution  of  the 
input  parameters  and  allows  the  simultaneous  identification  of  both  modeling  parameters 
and  MCR  properties  that  have  the  strongest  influence  on  the  peak  HGL  temperature 
achieved  during  our  scenario  [10]. 

The  results  of  our  case  study  indicate  that  for  existing  plants  the  one  controllable 
factor  available  to  an  operator  to  mitigate  the  probability  of  a  MCR  abandonment 
scenario  is  the  MCR  ventilation  rate. 

For  plants  yet  to  be  constructed,  the  effects  of  a  MCR  casualty  could  be 
substantially  reduced  through  an  improved  design  of  the  Remote  Shutdown  Panel  (RSP) 
that  provides  the  operators  with  the  necessary  cues  and  independent  circuitry  to  take 
actions  to  terminate  this  casualty  after  they  are  forced  to  abandon  the  MCR. 

The  results  of  this  study  are  heavily  dependent  on  the  distribution  of  the  Heat 
Release  Rate  (HRR)  of  the  prescribed  fire.  The  research  into  HRR  distributions  that  is 
currently  being  conducted  will  help  to  provide  regulators  with  the  necessary  tools  to  fully 
assess  the  contribution  to  core  damage  frequency  from  fires  initiated  from  within  the 
MCR. 
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II.  SCENARIO  IDENTIFICATION 


In  our  efforts  to  investigate  uneertainty  and  sensitivity  analysis  methods  for  risk 
scenarios  involving  binary  variables  and  mechanistic  codes,  we  have  chosen  to  use  a  fire- 
induced  accident  scenario  as  a  case  study. 

To  identify  a  fire-induced  risk  scenario  that  would  yield  interesting  results,  we 
considered  the  following  sources  of  uncertainty  specific  to  a  fire  PRA  [16]; 

1 .  Fire  Ignition  Frequency 

2.  Fire  Growth  and  Propagation 

3.  Fire  Suppression  Probability 

4.  Human  Error  Following  the  Fire  Event 

5 .  Mitigating  System  Availability 

Our  analysis  focused  on  fire-induced  scenarios  specific  to  a  plant  analyzed  as  part 
of  NUREG-1 150,  Severe  Accident  Risks:  An  Assessment  for  Five  U.S.  Nuclear  Power 

o 

Plants  [4],  which  identified  five  scenarios  with  a  core  damage  frequency  greater  than  10' 
yr'\  Eour  of  these  scenarios  contribute  to  over  99%  of  the  risk  of  core  damage  due  to 
fires  [17]. 

Three  of  the  scenarios  (fires  in  the  Emergency  Switchgear  Room  (ESWGR), 
Auxiliary  Building  (AUX  BEDG),  and  Cable  Vault/Tunnel  (CV/T))  are  similar  in  that  a 
fire  causes  the  loss  of  both  the  High  Pressure  Injection  (HPI)  and  Component  Cooling 
Water  (CCW)  systems.  The  loss  of  these  systems  leads  to  a  reactor  coolant  pump  (RCP) 
seal  Eoss  of  Coolant  Accident  (EOCA).  The  sequence  of  events  for  these  scenarios  is 
depicted  in  figure  2-1. 
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TURB 

TRIP 

W/MFW 


SCRAM 


-0 


NO  FAILURE 
OF  PORV 
TO  CLOSE 


NO  FAILURE 
OF  MFW/AFW 


■0 


HPI  SYS  FAIL 
IN  SEAL  MODE 


CCW  FAIL 
TO  COOL  RCP 


o 


HPI  SYS  FAILg 
- 0 


CORE  DAMAGE  ^ 
- <1 


Figure  2-1:  RCP  Seal  LOCA  Sequence  of  Events 


For  the  ESWGR  fire,  the  eore  damage  frequency  equation  is  given  in  NUREG- 
1150  as  follows: 

CDF  = 

Where: 

CDF  =  The  fire  induced  core  damage  frequency  for  the  ESWGR. 

^swgr  ^  The  frequency  of  ESWGR  fires  (of  all  sizes  and  severity). 

Qi'^G )  “  The  percentage  of  fires  in  the  data  base  where  the  fire  was  not 
manually  extinguished  before  the  COMPBRN  predicted  time  to 
critical  damage  occurred. 

=  The  probability  that  operators  will  fail  to  cross  connect  the  HPI  system 

prior  to  a  RCP  seal  LOCA.  This  does  not  require  operators  to  take 
action  in  the  direct  vicinity  of  the  fire. 

=  The  area  ratio  within  the  ESWGR  for  a  small  fire  where  critical 

damage  occurred.  This  is  calculated  by  dividing  the  area  in  the 
ESWGR  where  a  small  fire  could  damage  both  the  HPI  and  CCW 
systems  by  the  total  area  of  the  ESWGR. 
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/jj  =  The  severity  ratio  of  small  fires  (based  on  generic  combustible  fuel 
loading). 

/^2  =  The  area  ratio  within  the  ESWGR  for  a  large  fire  where  critical 

damage  occurred.  This  is  calculated  by  dividing  the  area  in  the 
ESWGR  where  a  large  fire  could  damage  both  the  HPI  and  CCW 
systems  by  the  total  area  of  the  ESWGR. 

/^2  =  The  severity  ratio  of  large  fires  (based  on  generic  combustible  fuel 

loading). 

Eor  the  AEIX  BEDG  fire,  the  core  damage  frequency  equation  is  given  in 
NUREG-1150  as  follows: 

Where: 

CDF  =  The  fire-induced  core  damage  frequency  for  the  AEIX  BEDG. 

\ux  ^  The  frequency  of  AUX  BEDG  fires. 

=  The  area  ratio  within  the  AUX  BEDG  where  critical  damage 
occurred. 

=  The  severity  ratio  for  a  large  fire  (based  on  generic  combustible  fuel 
loading). 

g(rg)  =  The  percentage  of  fires  within  the  suppression  data  base  where  the  fire 
was  not  manually  extinguished  before  the  COMPBRN  predicted  time 
to  critical  damage  occurred. 
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=  The  probability  that  operators  will  fail  to  cross  connect  the  HPI  system 

prior  to  a  RCP  seal  LOCA.  This  requires  the  operator  to  take 
action  in  the  direct  vicinity  of  the  fire.  Because  of  this,  no  recovery 
action  was  allowed  until  15  minutes  after  the  fire  was  extinguished. 

For  the  CV/T  fire,  the  core  damage  frequency  equation  is  given  in  NUREG-1 150 
as  follows: 

CDF=A^^JJMTG)Qau,oRo, 

Where: 

CDF  =  The  fire-induced  core  damage  frequency  for  the  CV/T. 

^csr  ^  The  frequency  of  CV/T  fires. 

=  The  area  ratio  within  the  CV/T  where  critical  damage 
occurred. 

=  The  severity  ratio  (based  on  generic  combustible  fuel  loading). 

Qi^G  )  “  The  percentage  of  fires  in  the  data  base  where  the  fire  was  not 
manually  extinguished  before  the  COMPBRN  predicted  time  to 
critical  damage  occurred. 

Qauto  ^  The  probability  of  the  automatic  CO2  system  not  suppressing  the  fire 
before  the  COMPBRN  predicted  time  to  critical  damage  occurred. 

R^p  =  The  probability  that  operators  will  fail  to  cross  connect  the  HPI  system 

prior  to  a  RCP  seal  LOCA.  This  does  not  require  operators  to  take 
action  in  the  direct  vicinity  of  the  fire. 
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The  next  scenario  is  a  fire  in  the  Charging  Pump  Service  Water  Pump  Room 
(CPSWPR)  that  is  caused  by  a  general  transient  followed  by  a  stuck  open  Power 
Operated  Relief  Valve  (PORV)  that  leads  to  a  small  LOCA.  The  sequence  of  events  for 
this  scenario  is  depicted  in  Figure  2-2. 


FAILURE  OF  HPI  SYS  FAIL^,  CORE  DAMAGE  ^ 

TURB  PORV  TO  CLOSE^J^  ^  ^ 

TRIP  SCRAM  9 

w/  MFW  r  'i  ' 

- CT  I - O 


Figure  2-2:  CPSWPR  Sequence  of  Events 


For  the  CPSWPR  fire,  the  core  damage  frequency  equation  is  as  follows: 

CDF  =  A^,.Q{rG)Qpor. 

Where: 

CDF  =  The  fire-induced  core  damage  frequency  for  the  CPSWPR. 

=  The  frequency  of  CPSWPR  fires. 

Qi'^G )  “  The  percentage  of  fires  in  the  data  base  where  the  fire  was  not 
manually  extinguished  before  the  COMPBRN  predicted  time  to 
critical  damage  occurred. 

Qporv  ~  The  probability  of  having  a  stuck  open  PORV  with  failure  to  isolate 
the  leak. 


The  final  fire  scenario  is  a  fire  in  benchboard  1-1  that  causes  a  spurious  Power 
Operated  Relief  Valve  (PORV)  lift  and  forced  Main  Control  Room  (MCR)  abandonment 
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followed  by  failure  to  reeover  the  plant  from  the  Remote  Shutdown  Panel  (RSP).  In  this 
scenario,  PORV  indication  is  not  provided  at  RSP  and  the  PORV  “disable”  function  on 
the  RSP  is  not  electrically  independent  from  the  MCR. 

A  fire-induced  manual  scram  or  turbine  trip  (external  event)  initiates  the  internal 
event  tree  T3-Q-D1.  The  sequence  begins  with  Figure  2-3. 

Where; 

T3;  Turbine  Trip  with  Main  Feed  Water  (MFW)  Available. 

Q:  Failure  of  a  PORV  to  close  after  transient  (GO  TO  S2). 

Dl;  Failure  of  charging  pump  system  in  High  Pressure  Injection  (HPI)  mode. 


RCI 

(Q) 


AFW 

(I-) 


MFW 

(M) 


SEAL 

COOL 

ccw 

HPI 

PRV 

CONT 

SYS 

CORE 
VULNR 
TO  CD 

LPR 

HPR 

(D3| 

.^(WL_ 

-0  T3 


-<  T3^)3 


-C  T3^)3-W 
-<|  T3-L 
-<  T3UM 


RPS 

(K) 


■6 


TURB 

TRIP 

W/MFP 


(T3) 


■6 


o - 0 - 6 


-6 


t31-mt;s 

T31-MT;S-H2 

T3l-M^S-H1 

T3UM^SW 

T3UMJ= 

T3UMT32 


<|  T3.Q 
C  T3T< 


Figure  2-3:  Event  Tree  T3-  Turbine  Trip  with  MFW  Available 


Once  Events  T3  and  Q  have  taken  place,  the  sequence  continues  on  Figure  2-4, 


the  event  tree  for  a  small  LOCA.  The  charging  system  does  not  initiate  in  HPI  mode 


because  the  stuck  open  PORV  does  not  cause  sufficient  depressurization. 


II 


OPER 
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(OP) 


LPR 

(HI) 


HPR 

(H2) 


HPI 

(D1) 


RPS 

(K) 


-6 


SMALL 

LOCA 

(S2) 
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-O- 
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(P1) 
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(CV) 


■o 


-6 


■o 


■o 


-o 


-6 


o 


{ 


■{ 


o- 

-o- 


i: 


-o 

-<i 

-<] 

-o 

-<3 

-<l 

-<3 

-<l 

-a 

-<l 

-<1 


-o 

-<3 

-<3 

-<] 

-O 


S2 

S2-H1 

S2-OD 

S2-OD-H2 

S2-OD-H1 

S2-CS 

S2-CS-H1 

S2-CS-CV 

S2-CS-OD 

S2-CS-OD-H2 

S2-CS-OD-H1 

S2-L 

S2-L-H2 

S2-L-H1 

S2-L-CS 

S2-L-CS-H2 

S2-L-CS-H1 

S2-L-CS-CV 

S2-L-P1 


-<3  S2-D1 


^  S2-K 


Figure  2-4:  Event  Tree  S2-Small  LOCA 

Figure  2-5  shows  how  the  fire  then  affects  individual  events  in  the  sequence: 


TURB 
TRIP 
W/  MFW 


SCRAM 


FAILURE  OF 
PORV  TO  CLOSE 


HPI  SYS  FAIL 


o- 


o 


o 


o 


CORE  DAMAGE  ^ 

- <1 


Figure  2-5:  Fire  Effects  on  Internal  Sequence  of  Events 
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For  the  MCR  fire,  the  eore  damage  frequeney  equation  given  in  NUREG-1 150  is 
as  follows: 

CDF  = 

Where: 

CDF  =  The  fire-indueed  eore  damage  frequeney  for  the  MCR. 

=  The  frequeney  of  MCR  fires. 

=  The  area  ratio  of  benehboard  1-1  to  total  eabinet  area  within  the  MCR. 
fj.  =  The  probability  that  operators  will  not  suoeessfully  extinguish  the  fire 
before  foreed  abandonment  of  the  MCR.  Aeeording  to  NUREG/CR- 
6850,  Fire  PRA  Methodology  for  Nuclear  Power  Facilities:  Volume  2: 
Detailed  Methodology  [2],  this  requires  either: 

(1)  Heat  flux  at  six  feet  above  the  floor  to  exeeed  1  kW/m  ,  based  on 
eausing  pain  to  skin,  whieh  equates  to  a  smoke  layer  of  95°C.  Or, 

(2)  The  smoke  layer  to  deseend  below  six  feet  from  the  floor  AND  the 
smoke  optieal  density  to  exeeed  3.0  m'\  eausing  the  operators  to  be 
unable  to  see  exit  signs. 

R^^  =  The  probability  that  operators  will  unsuoeessfully  reeover  the  plant 

from  the  RSP.  To  successfully  recover  the  plant,  the  operator  must 
shut  the  PORV  block  valve,  despite  not  having  indication  on  the  RSP 
that  the  PORV  has  lifted. 

Each  of  the  scenarios  under  consideration  would  provide  interesting  results  to  the 
fire  ignition  and  fire  growth  and  propagation  portions  of  our  analysis  due  to  the  increased 
data  and  computational  methods  available  to  update  the  CDE  distributions. 
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Evaluating  the  non-suppression  of  fire  events  is  typieally  done  through  statistieai 
methods  (Seetion  VIII)  involving  assessment  of  the  time  available  for  fire  suppression 
before  equipment  damage  or  an  operator  evaeuation  threshold  is  met.  These  statistieai 
models  are  often  based  on  time  eonstants  derived  from  historieal  fire  brigade 
performanee  data  and  automatie  system  reliability  data.  Beeause  the  system  reliability 
data  for  automatie  CO2  systems  had  to  be  modified  to  aeeount  for  the  short  time  to 
eritieal  damage  predieted  by  COMPBRN,  the  CV/T  seenario  may  not  yield  interesting 
results. 

The  CPSWPR  seenario  requires  a  signal  unrelated  to  the  fire  to  be  sent  to  lift  a 
PORV  and  the  subsequent  failure  of  the  PORV  to  reelose  and  isolate  the  leak.  This 
faetor  reduces  the  CPSWPR  scenario  Core  Damage  Frequency  (CDF)  contribution  to  less 
than  one  percent  of  the  overall  CDF  due  to  fires.  Because  of  this,  we  have  eliminated  the 
CPSWPR  scenario  from  further  consideration. 

The  ESWGR,  CV/T  and  CPSWPR  scenarios  do  not  require  human  failure  events 
under  increased  stress.  Only  the  AUX  BFDG  and  MCR  scenarios  rely  on  operator  action 
in  the  vicinity  of  the  fire  to  mitigate  the  sequence  of  events  leading  to  core  damage. 

The  ESWGR,  AUX  BFDG,  and  CV/T  scenarios  analyzed  in  NUREG-1 150  lead 
to  a  RCP  seal  FOCA.  Since  NUREG-1 150  was  published,  there  have  been 
improvements  implemented  in  RCP  seal  technology  [18,  19]  that  may  lower  the 
significance  of  re-analyzing  these  scenarios. 

Other  factors  we  considered  are  that  NUREG-1 824,  Verification  and  Validation 
of  Selected  Fire  Models  for  Nuclear  Power  Plant  Applications  [5],  specifically  covers  a 
MCR  fire  scenario  and  has  generated  experimental  data  relating  to  a  MCR  of  general 
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dimensions.  With  experimental  data  available,  a  Bayesian  methodology  for  determining 
model  uneertainty  appears  possible  [7]. 

As  we  move  forward,  we  intend  to  investigate  methods  to  analyze  uneertainty  and 
sensitivity  in  risk  seenarios  indueed  by  main  eontrol  room  fires. 
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III.  FIRE  IGNITION  FREQUENCY 


Since  NUREG-1 150  was  published  in  1990,  there  has  been  an  overall  downward 
trend  in  fire  ignition  frequency,  including  fires  initiated  in  the  MCR.  This  is  expected  as 
plants  have  improved  hre  prevention  policies  from  lessons  learned  and  knowledge 
sharing  practices.  Another  signifieant  factor  is  the  decline  in  cigarette  smoking  rates 
nationwide. 

We  will  now  develop  a  distribution  of  fire  ignition  frequeneies  speeific  to  our 
MCR  abandonment  scenario  by  applying  the  most  updated  data  and  methods  available. 

The  MCR  fire  ignition  frequeney  given  in  the  NUREG-1 150  analysis  is  a  Gamma 
distribution  characterized  by  a  shape  factor  a  =  1  and  scale  factor  P  =  555.56,  such  that 
the  probability  density  funetion  (pdf)  is  given  by  [20]; 


fix  \  a,p) 


/?“r(a) 


Where  r(a)  is  the  Gamma  function.  The  mean  value  of  the  Gamma  distribution  is 

—  1150  _3  _i 

aP'  .  Therefore,  the  mean  value  of  MCR  fire  ignition  frequency  is:  /I  cr  —  1.8x10' 

The  resulting  distribution  is  shown  in  Eigure  3-1. 


16 


MCR  Fire  Ignition  Frequency  Given  in  NUREG-1150 


O.OOE+OO  2.50E-03  5.00E-03  7.50E-03  1.00E-02 

yr*-1 

Figure  3-1:  MCR  Fire  Ignition  Frequency  Given  in  NUREG-1150 

EPRI  \Q\6135,  Fire  PRA  Methods  Enhancements:  Additions,  Clarifications,  and 
Refinements  to  EPRI  101189  [11]  was  published  in  2008  and  contains  the  latest  re- 
evaluation  of  fire  ignition  frequency  trends  and,  like  NUREG-1 150,  gives  the  MCR  fire 
ignition  frequency  as  a  Gamma  distribution  with  a  shape  factor  of  a  =1  but  differs  in  that 
a  scale  factor  of  P  =  1212.9  is  given.  This  revises  the  mean  predicted  value  of  MCR  fire 

—  EPRI  , 

ignition  frequency  to:  /I cr  =8.24x10  yr  and  results  in  the  distribution  given  in  Eigure 
3-2. 
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Updated  Generic  MCR  Fire  Ignition  Frequency 


O.OOE+OO  2.50E-03  5.00E-03  7.50E-03  1.00E-02 

yr^-l 

Figure  3-2:  MCR  Fire  Ignition  Frequency  Given  in  EPRI 1016735 

Using  updated  methods,  the  predicted  mean  MCR  fire  ignition  frequency  is  a 
factor  of  2.18  lower  than  previously  thought. 

EPRI  1016735  also  contains  guidance  on  frequency  estimation  parameters  that 
can  be  used  to  update  the  generic  fire  ignition  frequencies  for  the  plant  of  concern.  The 
most  significant  source  of  plant-to-plant  variability  appears  to  be  the  differences  in  event 
recording  and  reporting  practices.  As  the  practice  of  Fire  Probabilistic  Risk  Assessment 
progresses,  it  will  be  necessary  for  plants  to  strictly  adhere  to  standardized  reporting 
criteria. 

Simply  knowing  the  frequency  of  fires  that  occur  in  the  MCR  is  not  sufficient  for 
our  analysis.  For  the  scenario  under  consideration,  the  fire  must  occur  in  benchboard  1-1 
to  cause  the  PORV  to  lift.  All  of  the  event  fire  data  to  date  indicate  that  the  only  source 
of  MCR  fires  is  electrical  cabinets.  In  order  to  obtain  the  frequency  of  MCR  fires  that 
would  initiate  our  scenario,  an  area  ratio  was  developed  by  measuring  the  area  of 
benchboard  1-1  and  dividing  it  by  the  total  MCR  electrical  cabinet  area. 
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The  area  ratio  of  benchboard  1-1  to  total  cabinet  area  within  the  MCR  given  in  the 
NUREG-1 150  analysis  is  a  Maximum  Entropy  distribution  characterized  by  a  lower 
bound  of  a  =  0.028,  an  upper  bound  of  b  =  0. 12,  and  a  mean  of  |J  =  0.084  such  that  the 
probability  density  function  (pdf)  is  given  by  [21]; 

Where  fiiP  ^  0)  satisfies: 

_  be^  -  1 


resulting  in  the  distribution  given  in  Figure  3-3. 


Figure  3-3:  Area  Ratio  of  Benchboard  1-1  to  total  MCR  Cabinet  Area 


For  this  case  study,  it  was  not  necessary  to  update  the  area  ratio  given  in  NUREG- 
1 150  as  the  measurements  obtained  from  the  MCR  are  assumed  to  have  remained  the 
same. 
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In  order  to  obtain  the  overall  distribution  of  fires  in  benchboard  1-1,  it  was 
necessary  to  sample  from  the  MCR  fire  ignition  frequency  and  area  ratio  distributions. 
The  resulting  combined /I *  /„  distribution  is  given  in  Figure  3-4. 


Frequency  of  Fires  in  Benchboard  1-1 


yr^-i 


Figure  3-4:  Frequency  of  Fires  in  Benchboard  1-1 

The  mean  value  of  the  distribution  of  the  Frequency  of  Fires  in  Benchboard  1-1 
is:  /ic.  */«  =6.91x10 

Now  that  a  distribution  of  frequencies  of  initiating  events  has  been  developed,  we 
will  attempt  to  determine  a  set  of  environmental  conditions  in  the  MCR  as  a  fire  develops 
in  order  to  determine  the  probability  of  operator  abandonment  given  a  fire  in  benchboard 
1-1. 
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IV.  FIRE  SCENARIO  MODELING 


NUREG-1934,  Nuclear  Power  Plant  Fire  Modeling  Application  Guide  (FMAG)  [3],  is 
currently  a  draft  for  public  comment  that  contains  guidance  on  the  use  of  the  five  fire  models 
analyzed  as  a  part  of  NUREG-1824,  Verification  and  Validation  of  Selected  Fire  Models  for 
Nuclear  Power  Plant  Applications  [5].  In  modeling  our  MCR  abandonment  seenario,  we 
will  attempt  to  apply  the  guidanee  eontained  in  NUREG-1934  andNUREG/CR-6850, 
EPRl/NRC-RES  Fire  PRA  Methodology  for  Nuclear  Power  Facilities:  Volume  2: 

Detailed  Methodology,  TASK  11  [2].  Eaeh  of  these  doeuments  eontains  speeifie  guidanee 
on  the  modeling  of  MCR  fires. 

Modeling  Objectives 

In  modeling  the  seleeted  seenario,  our  objeetives  are  to  develop  an  evolving  set  of 
environmental  eonditions  in  the  MCR  after  the  start  of  a  fire  in  benehboard  1-1  in  order 
to: 

1 .  Assess  the  length  of  time  the  MCR  remains  habitable  in  the  absenee  of 
suppression  efforts  by  eomparing  environmental  eonditions  to  the  MCR  abandonment 
eriteria  [2]: 

•  Heat  flux  at  six  feet  above  the  floor  to  exeeed  1  kW/m  ,  based  on 

eausing  pain  to  skin,  whieh  equates  to  a  smoke  layer  of  95°C. 

OR 

•  The  smoke  layer  to  deseend  below  six  feet  from  the  floor  AND  the 

smoke  optieal  density  to  exeeed  3.0  m‘\  eausing  the  operators  to 

be  unable  to  see  exit  signs. 
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2.  Provide  input  to  detection  and  suppression  models  (Section  VIII).  In  cases 
where  a  MCR  abandonment  condition  is  predicted  to  occur  by  the  model  in  the  absence 
of  suppression  efforts,  a  probability  of  non-suppression  must  be  determined  through 
assessment  of  the  time  available  for  suppression  between  fire  detection  and  forced 
abandonment.  In  this  scenario  there  are  no  installed  lire  suppression  systems,  but  the 
MCR  is  continuously  manned.  However,  because  the  lire  takes  place  within  a 
benchboard,  no  credit  is  taken  for  prompt  detection.  Redundant  heat  detectors  are  located 
directly  above  the  fire  and  automatic  detection  is  assumed  to  occur  with  a  negligible 
failure  probability. 

3.  Provide  input  to  Human  Reliability  Analysis  (HRA)  models  (Section  IX).  If 
the  operators  fail  to  suppress  the  fire  and  are  forced  to  abandon  the  MCR,  they  will  be 
required  to  take  action  from  the  Remote  Shutdown  Panel  (RSP)  under  increased  stress 
and  with  fewer  available  indications  to  prevent  core  damage. 

Fire  Model  Selection 

To  identify  the  optimum  model  to  analyze  the  scenario  under  consideration,  the 
live  lire  models  verified  and  validated  in  NUREG-1824  are  considered; 

(1)  Fire  Dynamics  Tools  (FDTs) 

(2)  Fire-Induced  Vulnerability  Evaluation,  Revision  1  (FIVE-Revl) 

(3)  Consolidated  Model  of  Fire  Growth  and  Smoke  Transport  (CFAST) 

(4)  MAGIC 

(5)  Fire  Dynamics  Simulator  (FDS) 
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There  are  three  general  types  of  fire  models.  FDTs  and  FIVE-Rev  1  are  libraries 
of  engineering  caleulations,  CFAST  and  MAGIC  are  zone  models,  and  FDS  is  a 
eomputational  fluid  dynamies  model. 

Although  FDTs  and  FIVE-Rev  I  provide  Hot  Gas  Eayer  (HGE)  temperature 
results,  they  are  not  suitable  for  this  scenario  because  they  do  not  provide  smoke 
concentration  or  heat  flux  data  to  completely  evaluate  each  of  the  MCR  abandonment 
criteria. 

EDS  results  have  been  shown  to  be  comparable  to  CEAST  and  MAGIC, 
especially  in  simple  configurations,  but  are  computationally  expensive.  The  single 
compartment  MCR  we  will  be  modeling  is  a  sufficiently  simple  configuration  that  a  zone 
model  will  provide  adequate  results.  Choosing  a  zone  model  allows  for  much  larger 
Monte  Carlo  samples  to  be  reasonably  achieved  in  determining  model  input  parameter 
uncertainties  (Section  VI)  and  sensitivity  analyses  (Section  X). 

In  NUREG-I824,  CEAST  and  MAGIC  received  identical  validations  for  the 
outputs  of  interest  in  a  MCR  abandonment  scenario.  CEAST  was  chosen  over  MAGIC 
because  it  is  more  accessible  and  is  well  supported  by  the  National  Institute  of  Standards 
and  Technology  (NIST). 

Fire  Model  Description 

According  to  NUREG-I824,  Verification  and  Validation  of  Selected  Fire  Models 
for  Nuclear  Power  Plant  Applications,  Volume  5:  Consolidated  Fire  and  Smoke 
Transport  Model  (CFAST)  [9]:  “CEAST  is  a  two-zone  fire  model  that  predicts  fire- 
induced  environmental  conditions  as  a  function  of  time.  In  order  to  numerically  solve 
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differential  equations,  CFAST  subdivides  eaeh  eompartment  into  two  zones  that  are 
assumed  to  be  uniform  in  temperature  and  speeies  eoneentration,” 

The  two  publications  distributed  by  NIST  relevant  to  CFAST  are  the  CFAST 
Technical  Reference  Guide  [22],  which  explains  the  assumptions  and  physics  of  the 
model,  and  the  CFAST  User’s  Guide  [6],  which  explains  how  to  implement  the  model. 
Fire  Model  Input 

In  1985,  Factory  Mutual  and  Sandia  National  Laboratories  (FMSNL)  conducted  a 
series  of  tests  to  provide  data  for  use  in  validating  computer  fire  environment  simulation 
models,  specifically  MCR  scenarios  [8]. 

One  of  these  tests  (FMSNL  21)  was  conducted  in  a  MCR  mock-up  with  a  lire 
simulated  in  a  benchboard,  similar  to  the  scenario  under  evaluation.  FMSNL  21  was 
conducted  with  a  peak  Heat  Release  Rate  (HRR)  of  470  kW,  a  value  that  is  estimated  to 
exceed  the  peak  HRR  of  greater  than  92%  of  benchboard  fires.  In  addition,  the  FMSNL 
21  test  was  conducted  at  a  relatively  low  ventilation  rate  of  one  room  change  per  hour. 
Sensitivity  to  these  particular  input  parameters  is  analyzed  in  Section  X. 

Model  input  was  chosen  to  mimic  this  test  so  that  the  output  data  could  be 
compared  to  experimental  results  in  the  development  of  model  uncertainty  estimates 
(Section  V).  Inputs  specific  to  the  FMSNL  21  test  are  given  in  Appendices  A,  B,  and  C. 

CFAST  allows  the  user  to  specify  the  following  input  parameters: 

(1)  Ambient  Conditions 

(2)  Compartment  dimensions 

(3)  Construction  materials  and  material  properties 

(4)  Dimensions  and  positions  of  flow  openings 
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(5)  Mechanical  ventilation  specifications 

(6)  Sprinkler  and  detector  specifications 

(7)  Target  specifications 

(8)  Fire  properties 

NUREG-1824  states  that  the  most  important  input  factor  is  the  user  specified 
HRR.  Although  fire  growth  has  been  observed  to  follow  a  t^  growth  curve  [23],  the  input 
to  CFAST  was  linear.  This  and  other  assumptions  and  simplifications  contribute  to 
differences  between  predicted  and  observed  results.  Figure  4-1  shows  how  FIRR  was 
specified. 


FMSNL  21  Heat  Release  Rate 


Figure  4-1:  FMSNL  21  Heat  Release  Rate 

Fire  Model  Generated  Conditions 

Each  of  the  MCR  abandonment  criteria  (optical  density,  heat  flux,  and  FIGL 
temperature)  are  conditionally  dependent  on  the  HGL  descending  below  six  feet  from  the 
floor  so  that  the  conditions  actually  affect  the  operator. 
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CFAST  subdivides  each  compartment  into  two  zones  of  varying  volume,  an  upper 
HGL  and  a  lower  layer,  with  a  changing  interface  height  defined  as  the  HGL  height. 

The  first  step  in  our  analysis  was  to  ensure  that  the  MCR  operators  would  be 
exposed  to  the  conditions  predicted  in  the  HGL  layer  of  our  two  zone  fire  model  by 
determining  the  HGL  height  as  a  function  of  time  after  fire  ignition. 

As  shown  in  Figure  4-2,  the  HGL  height  descends  below  six  feet  in  just  under  five 
minutes. 


CFAST  HGL  Height 


Figure  4-2:  CFAST  HGL  Height 

Once  it  was  determined  that  the  HGL  would  in  fact  descend  to  a  level  low  enough 
to  affect  the  MCR  operators,  each  of  the  three  MCR  abandonment  criteria  were  evaluated 
as  a  function  of  time. 

Optical  density  is  a  measure  of  the  transparency  of  smoke.  It  depends  on  the  yield 
of  different  species  in  the  soot.  The  higher  the  optical  density,  the  lower  the  visibility: 

I/C 
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Where: 


I/Io  =  The  fraction  of  light  not  scattered  or  absorbed. 

X  =  The  optical  density  (units  of  length'^). 

X  =  The  straight  line  path  of  length  x  (units  of  length). 

In  order  for  the  optical  density  to  force  an  abandonment  condition,  the  HGL  must 
descend  below  six  feet  from  the  floor  (Figure  4-2)  and  the  smoke  optical  density  must 
exceed  3.0  m'\  While  CFAST  predicts  that  the  HGL  will  descend  from  the  ceiling 
relatively  quickly,  it  does  not  predict  that  the  optical  density  threshold  of  3.0  m''  will  be 
exceeded.  This  is  shown  in  Figure  4-3. 

CFAST  Optical  Density 
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Figure  4-3:  CFAST  Optical  Density 

In  addition  to  not  predicting  abandonment  due  to  optical  density,  CFAST  also 
does  not  predict  that  the  heat  flux  will  exceed  the  1000  W/m  threshold  to  force  the 
operators  to  abandon  the  MCR.  This  is  illustrated  in  Figure  4-4. 
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CFAST  Heat  Flux  to  Operator 
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Figure  4-4:  CFAST  Heat  Flux  to  Operator 

CFAST  predicts  that  the  MCR  will  become  uninhabitable  in  just  less  than  15 
minutes  due  to  the  HGL  temperature  and  height,  as  shown  in  Figure  4-5. 


CFAST  HGL  Temperature 
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Figure  4-5:  CFAST  HGL  Temperature 
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Modelins  Conclusions 


Even  at  the  relatively  low  ventilation  rate  prescribed  in  this  model,  the  optical 
density  is  only  predicted  to  reach  about  10%  of  the  necessary  value  to  force  evacuation. 
Although  the  heat  flux  from  the  smoke  layer  can  be  approximated  by: 

"  '"p  4 

(ir 


Where: 

=  The  radiated  heat  flux  (W/m^). 

(7  =  The  Stefan-Boltzmann  constant  (5.67x10'*  W/m^K"^). 

=  The  temperature  of  the  smoke  layer  (K). 

NUREG-1824  determined  that  the  CFAST  fire  model  is  more  capable  of 
accurately  predicting  HGL  layer  temperature  and  HGL  height  than  heat  flux. 

HGL  layer  height  is  predicted  to  descend  rapidly  after  fire  ignition,  while  the 
HGL  temperature  is  predicted  to  take  several  additional  minutes  to  reach  a  value  that 
would  force  evacuation. 

For  these  reasons,  I  intend  to  focus  on  the  HGL  temperature  criterion  (95°C)  for 
forced  MCR  abandonment. 


29 


V.  MODEL  UNCERTAINTY 


Fire  growth  and  propagation  contains  two  primary  sources  of  uncertainty.  The 
first  comes  from  the  input,  or  parameter,  uncertainty  that  occurs  due  to  the  distribution  of 
the  input  parameter  of  interest,  such  as  Heat  Release  Rate  (HRR)  (discussed  in  Section 
VI).  The  other  source  of  uncertainty  is  the  model  uncertainty,  which  is  the  epistemic 
uncertainty  resulting  from  the  assumptions  and  approximations  used  within  the  model. 

In  our  model,  CFAST,  the  primary  simplification  is  that  each  compartment  is 
divided  into  two  zones  that  are  assumed  to  have  uniform  properties.  Another 
simplification  is  made  by  not  solving  the  momentum  equation  explicitly.  However,  the 
conservation  of  mass  and  energy  equations  are  solved  as  ordinary  differential  equations. 

For  our  scenario,  the  primary  objective  of  assessing  model  uncertainty  is  to 
determine  the  probability  of  exceeding  the  MCR  abandonment  criterion  of  interest,  HGL 
temperature  greater  than  95  °C,  given  a  model  prediction. 

The  general  method  of  evaluating  model  uncertainty  is  through  the  comparison  of 
model  data  with  that  of  actual  experiments. 

A  Bayesian  Framework  for  Model  Uncertainty  Considerations  in  Fire  Simulation 
Codes  [7],  proposed  by  the  University  of  Maryland  (UMD),  was  first  conducted  by 
comparing  CFAST  output  to  experimental  results  contained  inNUREG-1824’s 
Benchmarking  Exercise  Three  (BE3)  [9].  Through  cooperation  with  the  United  States 
Nuclear  Regulatory  Commission  (NRC)  and  the  Massachusetts  Institute  of  Technology 
(MIT),  the  UMD  method  was  extended  to  include  data  from  the  EMSNL  21  and  22  tests 
conducted  as  part  of  NUREG/CR-4527,  An  Experimental  Investigation  of  Internally 
Ignited  Fires  in  Nuclear  Power  Plant  Control  Cabinets  Part  II:  Room  Effects  Tests  [8]. 
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The  results  from  the  UMD  method  using  only  the  BE3  data  (UMD  BE3)  are 
compared  to  the  updated  results  determined  by  using  the  EIMD  method  and  including  the 
EMSNE  21  and  22  test  data  (UMD  BE3  +  EMSNE  21/22).  The  results  from  the  UMD 
method  are  then  compared  to  the  method  presented  in  NUREG-1934  (EMAG)  [3]. 

Comparison  of  test  data  from  [8,  9]  to  model  predictions  show  that  CEAST 
consistently  over-predicts  HGE  temperature.  Eor  our  scenario,  the  discrepancy  is 
illustrated  in  Eigure  5-1. 


HGL  Temperature 
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Figure  5-1:  Comparison  of  Model  Prediction  to  Experimentally  Determined  HGL  Temperature 

The  assumptions  and  simplifications  made  by  our  model  lead  to  a  difference  in 
the  predicted  versus  experimentally  determined  values  that  is  especially  significant  in  the 
EMSNE  21  test  case  because  the  model  predicts  an  abandonment  condition,  while  the 
experiment  does  not. 

Another  factor  that  must  be  considered  is  the  uncertainty  present  in  the 
measurement  of  the  EMSNE  test  data.  The  model  uncertainty  techniques  presented  here 
take  this  uncertainty,  referred  to  as  experimental  uncertainty,  into  account. 
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Bayes’  theorem  allows  a  prior  probability  distribution  to  be  updated  to  posterior 
probability  distribution  when  new  evidenee  beeomes  available.  The  new  evidenee,  or 
data,  is  represented  by  a  likelihood  funetion,  sueh  that: 

L{data\  6)f^{6) 

J  \^)  r 

jL(data  \  6)f^{6)d6 

0 

Where: 

/ {6)  :  Posterior  probability  distribution. 

{6) :  Prior  probability  distribution. 

L{data  \  6) :  Likelihood  of  the  evidenee. 


The  UMD  method  [7]  uses  Bayesian  inferenee  to  update  a  prior  probability 
distribution  using  a  likelihood  funetion  derived  from  the  eomparison  of  fire  model  output 
to  experimental  data  as  follows: 
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Assuming  model  and  experimental  errors  are  independent  and  log-normally 
distributed,  the  likelihood  funetion  to  derive  the  posterior  joint  distribution  of  bmand  Sm 
beeomes: 

F  ^LN\b  -b  ,Js^  +s^  ] 

em  \  m  e  ’  \  m  e  ) 

Where: 

X:  Real  quantity  of  interest. 

Xm:  Model  predietion. 

Xei  Result  of  experiment. 

Fm!  Multiplieative  error  of  model  to  real  value. 

Fei  Multiplieative  error  of  experiment  to  the  real  value. 

Fem:  Multiplieative  error  of  experiment  to  model  predietion. 

bm:  Mean,  error  of  model  to  the  real  value. 

be:  Mean,  error  of  experiment  to  the  real  value. 

Sg:  Standard  deviation,  error  of  experiment  to  the  real  value. 

Sm:  Standard  deviation,  error  of  model  to  the  real  value. 

Onee  the  likelihood  funetion  has  been  developed,  the  posterior  joint  distribution 
of  bm  and  Smis  developed  as  follows: 

j  |/„(4,.0*AA  I  K,s„ 
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Where: 


L{X^,X^,b^,s^\b^,sJ  =  ^- 
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Where: 

fo  ’^m)-  Prior  joint  distribution  of  parameters. 

fm(bm,s„  I  X^i,X^.,b^,sJ  :  Posterior  joint  distribution  of  parameters. 


From  this,  a  distribution  of  the  real  quantity  of  interest  (X)  given  a  model 
prediction  (X^)  can  be  created  using  the  WinBUGS  (Microsoft  Windows  Bayesian 
Inference  Using  Gibbs  Sampling)  scripts  included  in  Appendices  D  and  E.  Such  that: 

F^-LN{b^,sJ 

X  = 

X-LN{\n{XJ  +  b^,sJ 


Figure  5-2  and  Table  5-1  show  the  results  of  the  UMD  method,  using  only  the 
BE3  data,  as  a  comparison  of  a  model  prediction  to  an  experimentally  determined  value 
with  both  experimental  and  model  uncertainty  boundaries  given. 
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The  scatter  of  data  is  assumed  to  result  from  uncertainty  in  both  the  model 
prediction  and  experimental  result.  Therefore,  neither  set  of  bounds  should  necessarily 
capture  the  entire  scatter  alone.  Also,  the  experimental  and  model  uncertainty  bounds  are 
for  the  real  value  given  a  model  prediction,  and  in  cases  like  ours  in  which  there  is  a  clear 
bias  in  the  model  prediction,  the  data  might  not  fall  within  the  bounds. 


CFAST  HGL  Temperature  Rise  BE3 


Figure  5-2:  CFAST  HGL  Temperature  Predictiou  vs.  BE3  Experimeutal  Data 
This  results  in  the  following  output  data; 


Fm=  0.899 


bm  =  -0.107 


Sm  =  0.02678 


Table  5-1:  UMD  BE3  WiuBUGS  Data 
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which  we  can  now  use  to  calculate  the  probability  of  exceeding  our  MCR  abandonment 
criterion  of  95  °C  given  a  model  prediction  (Xm): 


Pr(//GZ  >  95  I  Xm)  =  1  -  — exp[-^] 

{Itt)  1 


Where; 

^__ln(95)-h„,-ln(XJ 

Figure  5-3  illustrates  how  a  CFAST  model  prediction  (X^)  relates  to  the 
probability  of  exceeding  our  MCR  abandonment  criterion  of  95  °C. 

UMD  Probability  of  Exceeding  HGL  Temp  of  95  C 
(BE3  Only) 
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Figure  5-3:  UMD  BE3  Probability  of  Exceeding  HGL  Temp  of  95  C 

For  this  case  study,  the  UMD  Framework  was  then  extended  to  include  a 
comparison  of  CFAST  output  and  the  corresponding  FMSNL  tests  specific  to  this  MCR 
fire  scenario,  ft  is  important  to  note  that  the  BE3  experiments  and  the  FMSNL  test  series 
were  conducted  in  different  geometries  and  with  different  prescribed  fires.  Given  the 
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assumptions  contained  within  the  CFAST  model  and  its  known  sensitivity  to  HRR,  the 
additional  data  gained  from  the  FMSNL  test  series  may  or  may  not  be  expected  to  refine 
the  uneertainty  bounds  caleulated  using  only  the  BE3  test  data. 

Figure  5-4  and  Table  5-2  show  the  results  of  the  UMD  method,  using  both  the 
BE3  and  FMSNE  21/22  data,  as  a  comparison  of  a  model  prediction  to  an  experimentally 
determined  value  with  both  experimental  and  model  uncertainty  boundaries  given. 


CFAST  HGL  Temperature  Rise  BE3  +  FMSNL21/22 


Experimental  Value  (°C) 


Figure  5-4:  CFAST  HGL  Temperature  Predictiou  vs.  BE3  aud  FMSNL  21/22  Experimeutal  Data 


This  results  in  the  following  output  data; 
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Fm=  0.873 


bm  = -0.1614 


Sm  =  0.217 


Table  5-2:  UMD  BE3  and  FMSNL  WinBUGS  Data 


which  we  can  now  use  to  calculate  the  probability  of  exceeding  our  MCR  abandonment 
criterion  of  95  °C  given  a  model  prediction  (Xm).  This  is  illustrated  in  Figure  5-5. 


UMD  Probability  of  Exceeding  HGL  T emp  of  95  C 
(BE3  +  FMSNL  21/22) 


Figure  5-5:  UMD  BE3  Probability  of  Exceeding  HGL  Temp  of  95  C 


Although  adding  the  two  data  points  from  the  FMSNL  21  and  22  tests  caused  the 
uncertainty  bounds  predicted  by  the  UMD  method  to  widen,  it  is  important  not  to 
compare  these  methods  based  on  the  probability  of  exceedance  of  a  single  model 
prediction  only.  In  Section  VII  we  will  compare  these  methods  by  coupling  Monte  Carlo 
simulations  to  the  deterministic  CFAST  model  to  determine  an  integrated  probability  of 
exceedance  of  the  MCR  abandonment  criterion  of  interest. 

The  FMAG  method  uses  the  results  of  NUREG-1824  to  give  a  bias  factor  and 
model  error  for  each  quantity  of  interest  and  fire  model  used  while  assuming  model  error 
is  normally  distributed  such  that; 
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X  =  N 


S 


s 


) 


Where; 


X:  Real  quantity  of  interest. 

Xml  Model  prediction. 

5:  Bias  factor. 

;  Relative  model  error. 

For  the  CFAST  model  predicting  HGL  temperature  the  FMAG  gives  the 
following  data: 


6  =  1.06 


cr„ 


0.12 


Table  5-3:  FMAG  Data  for  CFAST  Predicting  HGL  Temperature 


For  the  FMSNL  21  test  case  modeled  in  CFAST  and  predicting  a  peak  HGL 
temperature  of  103  °C,  this  results  in  a  normal  distribution  with  a  mean  of  98.02  °C  and  a 
standard  deviation  of  9.96  °C  as  shown  in  Figure  5-6. 
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FMAG  Method:  Distribution  of  True  HGL  Temp  given  a  103  C 

CFAST  Prediction 


0  25  50  75  100  125  150  175  200 

Temperature  (C) 

Figure  5-6:  FMAG  Distribution  of  True  HGL  Temp 

We  can  now  use  this  data  to  calculate  the  probability  of  exceeding  our  MCR 
abandonment  criterion  of  95  °C  given  a  model  prediction  (X^)  and  an  ambient 
temperature  of  15 

?x{HGL  >  95  I  Xm)  =  1  -  — ^^exp[-^] 

{Itv)  l 


Where: 


95-// 

cr 


//  =  15  + 


(Xm-15) 

d 


cr  =  cr^  =*= 


^  Xm-\S^ 
~d 


Figure  5-7  illustrates  how  a  CFAST  model  prediction  (Xm)  relates  to  the 
probability  of  exceeding  our  MCR  abandonment  criterion  of  95  °C. 
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Figure  5-7:  FMAG  Probability  of  Exceeding  HGL  Temp  of  95  C 


From  these  results,  it  appears  as  though  the  method  presented  in  the  FMAG 
produces  results  that  are  more  conservative  than  the  UMD  method.  However,  at  this 
stage  it  would  be  premature  to  select  one  of  the  above  methods  of  determining  model 
uncertainty.  As  we  move  forward  we  will  compare  how  each  of  the  available  methods 
affects  the  outcome  of  our  scenario  from  a  probabilistic  standpoint  by  coupling  Monte 
Carlo  simulations  to  our  deterministic  fire  model,  CFAST,  to  account  for  the 
uncertainties  present  due  to  the  variance  in  the  distribution  of  input  parameters. 


41 


VI.  PARAMETER  UNCERTAINTY 


The  uncertainty  of  a  model  prediction  is  not  only  dependent  on  the  assumptions 
and  simplifications  of  the  model  itself,  but  also  on  how  the  uncertainties  in  input 
parameters  are  propagated  throughout  the  model.  Inputs  to  deterministic  fire  models  are 
often  not  precise  values,  but  instead  follow  statistical  distributions. 

In  simple  cases,  when  the  effect  on  a  single  output  quantity  due  to  changing  just 
one  input  parameter  is  desired,  empirical  correlations  may  be  appropriate.  NUREG-1934 
[3]  offers  very  useful  guidance  on  the  use  of  model-independent  empirical  correlations 
that  provide  one-to-one  mapping  of  the  effect  on  a  specific  output  quantity  given  a 
change  in  a  single  input  parameter.  For  our  scenario,  where  HRR  is  the  most  significant 
input  parameter  [9]  and  the  HGL  temperature  is  the  output  quantity  we  desire  to 
calculate,  the  McCaffrey,  Quintiere,  and  Harkleroad  (MQH)  correlation  [24]  is  given; 

Thgl-Thgl^  =CHHRRf'' 


Where; 

T  = 

Hot  Gas  Layer  Temperature. 

T  = 

tiGL„ 

Initial  Hot  Gas  Layer  Temperature. 

c 

Constant. 

HRR  = 

Heat  Release  Rate. 

The  value  of  the  constant  is  irrelevant,  as  the  relationship  we  are  seeking  is 
developed  by  differentiating  the  MQH  correlation  with  respect  to  HRR; 

_  2  MiRR 

'^BGL  ~'^HGL„  3  HRR 
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Where; 


^T, 


HGL 


T  -T 

^  HGL  ^  HGL„ 


Relative  Change  in  HGL  Temperature  Output. 


-  =  Relative  Change  in  HRR  Parameter  Input. 

HRR 

From  the  MQH  eorrelation,  it  is  expected  that  a  15%  increase  in  HRR  would  lead 
to  a  10%  increase  in  HGL  Temperature.  Figure  6-1  shows  a  comparison  the  FMSNL  21 
CFAST  simulation  with  HRR  increased  by  15%  compared  to  the  value  predicted  by  the 
MQH  correlation. 


Model  Prediction  vs.  Correlation 


Figure  6-1:  Comparison  of  CFAST  Output  to  MQH  Prediction 


From  Figure  6-1  it  is  apparent  that  the  MQH  correlation  is  useful  when  only  one 
input  parameter  is  subject  to  variability. 
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Due  to  the  eomplexity  and  non-linear  nature  of  our  fire  model,  empirieal  methods 
to  estimate  uneertainty  propagation  do  not  yield  suffieiently  refined  results  when  multiple 
input  parameters  are  allowed  to  vary.  In  order  to  more  adequately  assess  the  distributions 
of  fire  model  output  variables,  Monte  Carlo  simulations  have  been  eoupled  to  our  fire 
model,  CFAST,  through  a  tool  ealled  Probabilistie  Fire  Simulator  (PFS)  [10].  PFS  also 
gives  the  sensitivity  of  the  output  variables  to  the  input  variables  in  terms  of  rank  order 
eorrelation  eoeffieients  (Seetion  X). 

We  are  primarily  eoneemed  with  determining  the  probability  that  the  peak  HGL 
temperature  reaehed  in  our  MCR  fire  seenario,  as  predieted  by  CFAST,  will  exeeed  the 
95  °C  threshold  for  MCR  abandonment  in  the  absenee  of  suppression  efforts.  A 
seeondary  goal  of  interest  is  to  assess  the  time  available  for  suppression  efforts  by 
eomparing  the  deteetor  aetivation  time  with  the  time  of  foreed  abandonment.  PFS 
provides  the  neeessary  time  series  data  for  this  evaluation. 

Onee  the  probability  that  an  abandonment  eondition  in  the  absenee  of  suppression 
efforts  has  been  determined,  the  parameter  uncertainty  and  model  uncertainty  results  will 
be  combined  (Section  VII)  to  determine  an  integrated  probability  of  exceeding  the 
specified  abandonment  criterion  in  the  absence  of  suppression  efforts. 

Suppression  efforts  will  then  be  factored  in  through  an  assessment  of  the 
evolution  of  environmental  conditions  in  the  MCR  with  time  compared  to  detector 
activation  data  (Section  VIII).  From  this,  an  estimate  of  operator  abandonment  of  the 
MCR,  given  a  fire  in  benchboard  I -I,  will  be  obtained. 

Prior  to  fully  implementing  PFS,  it  was  necessary  to  first  use  the  spreadsheet 
environment  to  duplicate  the  modeling  results  obtained  in  Section  IV  in  order  to 
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demonstrate  proper  operation.  The  PFS  input  speeifieally  designed  to  mimie  the  FMSNL 
21  test  is  given  in  Appendix  F. 

Onee  fully  implemented,  PFS  eouples  Monte  Carlo  simulations  with  CFAST  in  a 
spreadsheet  environment  that  allows  the  user  to  vary  input  parameters  with  individually 
speeified  distributions. 

From  NUREG- 1824  [9]  it  is  known  that  CFAST  is  espeeially  sensitive  to 
variations  in  HRR.  In  order  to  adequately  model  the  FIRR  distribution  during  a  fire,  both 
the  initial  growth  period  and  the  fully  developed  state  (HRRmax)  must  be  eonsidered. 

Fire  growth  has  been  observed  to  follow  a  t^  growth  eurve  [23]  up  to  a  fully 
developed  state  where  the  HRR  beeomes  equal  to  HRRmax,  sueh  that: 


HRR{t)  =  min  j  4  000  *  — 


Where: 

HRR(t)  =  Heat  Release  Rate  as  a  funetion  of  time. 

HRRmax  =  Maximum  HRR  attained  during  the  fire  seenario. 
t  =  Time  after  fire  ignition, 

tg  =  HRR  growth  time. 

For  our  seenario,  tg  is  treated  as  a  normally  distributed  random  variable  with  a 
mean  of  320  seeonds  and  a  standard  deviation  of  100  seeonds,  eonsistent  with  [10]. 
Figure  6-2  shows  how  variations  in  tg  affeet  the  HRR  distribution  of  the  FMSNL  21  ease. 
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HRR  Growth  Time  Comparison 


Figure  6-2:  Comparison  of  HRR  Distributions  Given  Variations  in  Growth  Time 


The  maximum  value  that  HRR  reaehes  in  eaeh  simulation  (HRRmax)  used  in  this 
ease  study  is  given  in  Appendix  G  of  NUREG-6850  [2]  as  a  Gamma  distribution 
eharaeterized  by  a  shape  faetor  a  =  0.7  and  seale  faetor  P  =  216,  sueh  that  the  probability 
density  funetion  (pdf)  is  given  by  [20]; 


/(x  \  a,f5) 


Where  r(a)  is  the  Gamma  funetion.  The  mean  value  of  the  distribution  is  a*p. 
The  mean  value  of  HRRmax  given  in  Appendix  G  of  NUREG-6850  is; 

- 6850 

HRR  =151 .2kW  resulting  in  the  distribution  given  in  Eigure  6-3. 
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Heat  Release  Rate  Given  in  NUREG-6850 


0  200  400  600  800  1000 

kW 


Figure  6-3:  NUREG-6850  Appendix  G  HRR  for  Cabinets  with  Qualified  Cable 


Once  the  HRR  distribution  input  to  PFS  was  speeified,  it  was  also  determined  that 
we  would  vary  the  environmental  and  thermodynamie  properties  that  eould  have  an 
outcome  on  the  fire  scenario  analysis.  Table  6-1  lists  the  random  variables  used  in  our 


seenario  with  their  assoeiated  distributions  speeified  as  input  to  PFS. 


VARIABLE 

DISTRIBUTION 

UNITS 

Heat  Release  Rate 

Gamma 

alpha  =  0.7 

beta  =  216 

kW 

HRR  Growth  Time 

Normal 

mean  =  320 

s.d.  =  100 

M 

Ventilation  Rate 

Uniform 

min  =  0.38 

max  =  1.90 

m^3/s 

Lower  Oxygen  Limit 

Uniform 

min  =  9.00 

max  =  1 1 .00 

% 

Relative  Humidity 

Uniform 

min  =  45.00 

max  =  55.00 

% 

Radiative  Fraetion 

Uniform 

min  =  0.315 

max  =  0.3850 

Heat  of  Combustion 

Uniform 

min  =  4.05e4 

max  =  4.95e4 

(kJ/kg) 

Wall  Conduetivity 

Uniform 

min  =  0.108 

max  =  0.132 

W/m*K 

Wall  Speeifie  Heat 

Uniform 

min  = 1250 

max  =  1375 

J/kg*K 

Wall  Density 

Uniform 

min  =  720 

max  =  2200 

kg/m^3 

Wall  Thiekness 

Uniform 

min  =  0.0225 

max  =  0.275 

M 

Wall  Emissivity 

Uniform 

min  =  0.855 

max  =  1.045 

Table  6-4:  Random  Variables  Used 
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For  each  simulation,  a  CFAST  input  file  is  created  from  both  the  fixed  data  and 
the  samples  drawn  from  random  variables  within  the  above  constraints.  CFAST  is  then 
run  to  generate  and  save  time  series  data  of  the  user’s  choice  for  each  simulation.  For  our 
scenario  it  was  necessary  to  generate  detector  activation  time  data,  FIGL  temperature  time 
series  data,  and  a  record  of  the  peak  FIGL  temperature  achieved  in  each  simulation. 

Convergence  criteria  were  set  so  that  the  sample  size  (n)  would  be  sufficiently 
large  to  obtain  meaningful  results  while  remaining  computationally  inexpensive.  For  our 
scenario,  we  were  interested  in  determining  the  probability  that  that  a  given  fire  would 
have  a  predicted  peak  HGL  temperature  in  excess  of  95  °C  within  the  following 
constraints; 

Pr(r  >  95Cy  -  Pr(T  >  ^ 

Pr(r  >  95C)”^‘°“° 


Reaching  our  convergence  criteria  required  10,000  samples,  as  illustrated  in 
Figure  6-4. 


0.0500 

0.0400 

0.0300 

Pr(l>95C) 

0.0200 

0.0100 

0.0000 


Figure  6-4:  Convergence  of  Probabilistic  Fire  Simulator  Results 


Convergence  of  PFS  Simulations 
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The  distribution  of  peak  HGL  temperature  predictions  by  PFS  is  shown  in  Figure 


6-5  as  a  probability  mass  function  (pmf). 


PFS  Predicted  Peak  HGL  Temperature  Distribution 


Model  Prediction  (C) 

Figure  6-5:  Peak  Predicted  HGL  Temperature  Distributiou  Predictiou 

Once  the  distribution  of  peak  HGL  temperatures  is  known,  the  probability  that  an 
abandonment  condition  will  be  reached  can  be  derived  from  the  cumulative  distribution 
function  (CDF)  of  peak  predicted  HGL  temperature  given  in  Figure  6-6,  such  that: 


CDF  —  J  p^f  (T 

0 
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The  probability  of  reaching  an  abandonment  condition  is  then  calculated  as: 

95 

Pr(r  >  95C)  =  1  -  j 

0 

This  results  in; 

Pr(r>  950  =  0.0395 

This  quantity  represents  a  conservative  estimate  of  the  probability  of  operator 
abandonment  of  the  MCR  given  a  fire  in  benchboard  1-1.  In  Section  V  it  was  determined 
that  our  CFAST  fire  model  consistently  over-predicted  peak  HGL  temperature.  In  order 
to  refine  this  estimate  it  will  be  necessary  to  combine  model  and  parameter  uncertainties 
(Section  VII)  and  account  for  suppression  efforts  (Section  VIII)  taken  by  the  operators  in 
the  MCR. 
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VII.  COMBINED  MODEL  AND  PARAMETER  UNCERTAINTY 


Now  that  we  have  developed  separate  teehniques  for  assessing  model  and 
parameter  uncertainties,  our  next  task  is  to  combine  the  methods  in  an  effort  to  calculate  a 
cumulative  probability  of  exceeding  our  MCR  abandonment  threshold  given  a  fire  in 
benchboard  1-1  and  assuming  no  suppression  efforts  are  made. 

If  the  only  input  parameter  that  we  intended  to  vary  was  HRR,  and  if  the  HRR 
distribution  was  known  and  normally  distributed,  we  could  combine  the  NUREG-1934 
(FMAG)  model  uncertainty  (also  normally  distributed)  and  parameter  uncertainty  via 
quadrature,  such  that  [3]; 

Where: 

a  =  The  standard  deviation  of  the  combined  error. 

=  The  standard  deviation  of  the  model  error. 

p  =  A  sensitivity  factor  (2/3  for  HGL  temperature). 

aj  =  The  standard  deviation  of  the  input  (parameter)  error. 

However,  the  FMAG  method  is  not  applicable  to  our  case  study  because  the 
parameter  that  we  are  most  sensitive  to,  HRR,  is  not  normally  distributed  nor  is  it  the 
only  input  parameter  to  be  varied. 

In  Section  V  we  presented  model  uncertainty  techniques  that  provide  a  probability 
of  exceeding  a  threshold  given  a  single  model  prediction,  Xm.  We  will  now  extend  this 
framework  to  handle  a  distribution  of  model  predictions  (Figure  6-5)  that  is  created  when 
input  parameters  are  allowed  to  vary  as  specified  in  Section  VI  (Table  6-1). 
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Each  of  the  model  uncertainty  techniques  covered  in  Section  V  suggest  that  the 
true  value  that  peak  HGL  temperature  reaches  is  less  than  that  predicted  by  CFAST 
(being  run  by  PFS).  Figure  7-1  shows  a  comparison  of  our  model  uneertainty  teehniques 
by  the  mean  value  of  the  peak  HGF  temperature  predicted  in  each  model  simulation. 


Comparison  of  Model  Uncertainty  Techniques  by  Mean  Peak  HGL 
Temperature  Prediction 


Figure  7-1:  Comparison  of  Model  Uncertainty  Techniques  by  Mean  Peak  HGL  Temperature 
Although  CFAST  generally  over-predicts  peak  HGL  temperature,  evaluating  the 
probability  that  operators  will  be  forced  to  abandon  the  MCR  from  the  mean  predictions 
alone  is  insufficient.  Limiting  our  analysis  in  this  way  neglects  the  variance  in  the  model 
uneertainty  techniques  that  allows  for  the  possibility  that  a  real  value,  X,  may  exceed  a 
model  prediction,  Xm. 

To  generate  a  cumulative  probability  of  exeeeding  our  abandonment  threshold,  it 
was  necessary  to  consider  the  sum  of  the  contributions  from  each  simulation  and 
normalize  the  result  for  our  sample  size,  such  that: 
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10000 

XPr(7’“>95C|;\r„,,) 

CDF=  Pr(r"";  >  95C)  =  - 

HGL 

Ordering  our  model  predictions  by  ascending  peak  HGL  temperature,  the  model 
uncertainty  techniques  are  compared  by  cumulative  probability  of  exceeding  the 
abandonment  threshold  in  Figure  7-2. 


Cumulative  Distribution  of  the  Probability  of  Forced  MCR  Abandonment 


9000  9100  9200  9300  9400  9500  9600  9700  9800  9900  10000 

simulations 

(ordered  by  ascending  peak  HGL  temperature  model  prediction) 

Figure  7-2:  CDF  of  Forced  MCR  Abandonment 
It  is  apparent  from  Figures  7-1  and  7-2  that  the  most  conservative  method  of 
analyzing  our  scenario  is  to  neglect  model  uncertainty  altogether  and  use  the  values 
predicted  by  CFAST  (PFS).  Less  conservative  results  can  be  obtained  by  using  the 
method  presented  in  the  FMAG,  with  the  UMD  method  yielding  the  least  conservative 
results. 

Although  re-analyzing  the  UMD  method  with  the  inclusion  of  the  FMSNL  test  21 
and  22  data  caused  a  greater  variance  in  the  expected  value  of  interest,  the  aggregate 
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effect  of  this  difference  is  nearly  negligible  when  a  cumulative  probability  of  exceedance 
is  calculated  (Figure  7-3). 


Comparison  of  Model  Uncertainty  Techniques  by  Cumulative 
Probability  of  Exceedance 


Pr(l>95C|no 

suppression) 


Figure  7-  3:  Comparison  of  Model  Uncertainty  Techniques  by  Probability  of  Exceedance 
Rather  than  constrain  our  analysis  to  one  of  the  above  techniques  as  we  move 
forward,  comparative  results  will  be  offered.  However,  because  of  the  close  agreement  in 
the  UMD  methods  their  results  will  be  considered  as  one. 

While  we  are  closer  to  assessing  the  probability  that  operators  will  be  forced  to 
abandon  the  MCR  and  take  actions  from  the  RSP  in  case  of  fire  in  benchboard  1-1,  our 
results  are  still  conservative  in  that  no  credit  has  been  taken  for  manual  fire  suppression 
efforts.  In  Section  VIII  we  will  take  operator  suppression  efforts  into  consideration  as  we 
refine  our  estimate  of  the  probability  of  forced  operator  abandonment  of  the  MCR. 
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VIII.  FIRE  SUPPRESSION  ANALYSIS 


In  order  for  the  seenario  we  are  eonsidering  to  take  plaee  a  fire  must  eause  the 
MCR  to  be  abandoned.  To  this  point  we  have  only  eonsidered  whether  a  fire  could  eause 
an  abandonment  eondition,  in  this  seetion  we  eonsider  whether  a  fire  will  eause  an 
abandonment  eondition. 

In  evaluating  the  eombined  model  and  parameter  uneertainties  present  in  our 
seenario,  the  goal  was  to  develop  an  expression  for  the  probability  that  a  fire  in 
benchboard  1-1  would  force  operator  abandonment  of  the  MCR  if  no  suppression  efforts 
were  made.  The  final  step  in  our  fire  growth  and  propagation  analysis  is  to  assess  the 
probability  that  operators  in  the  MCR  will  suppress  a  fire  before  an  abandonment 
condition  occurs.  The  probability  of  forced  MCR  abandonment  is: 

>  95C  I  Non  -  Suppression)  *  Vv(Non  -  Suppression) 

To  ascertain  the  probability  of  non-suppression  an  event  tree  analysis  will  be 
conducted  using  the  method  presented  in  NUREG-6850  [2]. 

The  MCR  is  a  continuously  occupied  space  with  two  heat  detectors  located  on  the 
ceiling  such  that  once  a  fire  is  detected  suppression  efforts  will  begin  immediately. 
However,  there  are  no  automatically  actuated  fixed  suppression  systems  in  the  MCR. 

In  developing  the  event  tree  for  our  scenario,  no  credit  was  taken  for  prompt 
manual  detection  of  the  fire  due  to  its  location  within  a  benchboard  and  the  presence  of  a 
forced  ventilation  system  that  each  contributes  to  likelihood  that  the  operators  will  not 
detect  the  fire.  It  was  also  assumed  that  delayed  manual  detection  was  unnecessary  due 
to  the  redundant  heat  detectors  that  will  provide  automatic  detection  with  a  negligible 
failure  rate.  The  suppression  event  tree  for  our  scenario  is  given  in  Figure  8-1. 
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Figure  8-1:  Fire  Suppression  Event  Tree 


Through  an  assessment  of  manual  fire  suppression  historieal  data,  generie 
industry-wide  response  rates  (k)  for  different  plant  locations  have  been  developed  to  aid 
in  the  evaluation  of  the  probability  that  a  fire  will  be  manually  suppressed  as  a  function  of 
time.  From  the  model  presented  in  [2],  the  probability  that  a  fire  will  not  be  suppressed  is 
given  as: 


Pr(iV5)  = 


*  Suppression 


Where: 

Fr{NS)  =  The  probability  of  non-suppression. 


Suppression  Abandonment  Detection 

^MCR  =0.33min  '  =0.0055^  ' 


Finding  the  time  available  for  suppression  is  done  through  analysis  of  the  time 
series  data  generated  in  Section  VI.  Figure  8-2  is  an  example  of  how  the  time  available 
for  detection  is  obtained  using  the  FMSNL  21  test  case  as  an  example. 
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Suppression  Time  for  FMSNL  21  Test  Case 


0  300  600  900  1200  1500  1800 

time  (s) 

Figure  8-2:  FMSNL  Time  Available  for  Suppression 
To  develop  an  aeeurate  estimate  of  the  time  available  for  suppression  in  our 
seenario,  each  of  the  Monte  Carlo  samples  performed  by  PFS  (Section  VI)  that  predicted 
a  peak  HGL  temperature  greater  than  the  abandonment  criterion  were  considered.  This 
suppression  time  data  was  used  to  generate  Figure  8-3,  which  evaluates  each  sample  for 
the  probability  of  non-suppression. 


Probability  of  Non-Suppression 


Figure  8-3:  Probability  of  Non-Suppression 
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Once  the  probability  of  non-suppression  of  each  scenario  is  known,  the  mean  is 


ealeulated  as: 

395 

_  ElTf®), 

Fr{NS)  =  - =  0.3231 

395 

We  ean  now  estimate  the  mean  probability  of  operator  abandonment  of  the  MCR 
due  to  HGL  temperature  from  a  fire  in  benchboard  1-1  for  each  of  our  model  uncertainty 
techniques  as: 

=  0.0127 

fFMAG^  0.0111 

fUMD  ^  QQQ9J 

This  is  a  substantially  lower  estimate  than  that  provided  in  NUREG-1 150  [4] 
which  gives  the  probability  of  operator  abandonment  of  the  MCR  as  Maximum  Entropy 
distribution  charaeterized  by  a  lower  bound  of  a  =  0.01,  an  upper  bound  of  b  =  0.25,  and 
a  mean  of  p  =  0.1  sueh  that  the  probability  density  funetion  (pdf)  is  given  by  [21]: 

f(e\a,h,M)=  J, 

Where  ^  0)  satisfies: 

_  be'^  -  ae'^  1 

p 

This  results  in  the  distribution  given  in  Figure  8-4. 
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Probability  of  Abandonment  Given  in  NUREG-1150 


Figure  8-4:  MCR  Abandonment  Probability  Given  in  NUREG-1150 

MCR  abandonment  probability  has  been  an  issue  of  great  eontention  in  the  Fire 
PRA  eommunity.  Applieation  of  state-of-the-art  methods  in  estimating  the  probability  of 
operator  abandonment  in  this  seenario  suggest  that  the  value  given  in  NUREG-1 150  is 
overly  conservative;  however,  the  value  is  still  significant  and  MCR  abandonment 
scenarios  should  not  be  automatically  eliminated  from  consideration  when  conducting  a 
fire  PRA. 

In  Section  X  we  will  consider  factors  within  an  operator’s  control  to  mitigate  the 
probability  of  a  MCR  abandonment  scenario  occurring. 

Once  operators  are  forced  to  abandon  the  MCR  they  are  required  to  take  actions 
from  the  RSP  to  prevent  core  damage.  The  probability  of  successful  operator  action  will 
be  considered  in  the  next  section. 
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IX,  HUMAN  RELIABILITY  ANALYSIS 


After  the  fire  in  benehboard  1-1  eauses  spurious  aetuation  of  a  PORV  and  the 
operators  are  foreed  to  abandon  the  MCR  and  station  themselves  at  the  RSP,  they  must 
shut  the  PORV  bloek  valve  to  prevent  eore  damage.  A  low  probability  of  sueeessful 
operator  aetion  is  expeeted  due  to  the  PORV  closure  status  not  being  displayed  on  the 
RSP  and  the  PORV  block  valve  controls  on  the  RSP  not  being  electrically  independent  of 
the  MCR  benchboard  where  the  fire  is  occurring. 

NUREG-1921,  EPRI RES  Fire  Human  Reliability  Analysis  Guidelines  [12],  is 
currently  a  draft  for  public  comment  that  provides  guidance  on  estimating  Human  Error 
Probabilities  (HEPs)  for  Human  Eailure  Events  (HEEs).  Three  general  methods  are  given 
for  analyzing  post-fire  human  error  probability: 

(1)  Screening  HRA  Quantification 

(2)  Scoping  HRA  Quantification 

(3)  Detailed  HRA  Quantification 

These  methods  are  sequentially  less  conservative  and  more  detailed. 

The  screening  method  is  the  simplest  and  most  conservative  of  the  methods  and  is 
often  used  as  a  first  step  in  determining  which  sequences  warrant  further  consideration. 
The  scoping  method  is  less  conservative  than  the  screening  method  and  uses  decision  tree 
logic  to  assign  appropriate  HEPs.  Detailed  HEP  quantification  can  be  accomplished 
through  the  use  of  the  EPRI  HRA  CAECUEATOR  [25]  or  NUREG-1880,  A  Technique  for 
Human  Event  Analysis  (ATHEANA)  [26].  Both  of  the  detailed  quantification  methods 
require  in-depth  analysis  of  plant  specific  training  information  and  procedures  that  are 
beyond  the  scope  of  this  generalized  report. 
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The  screening  method  of  quantifying  HEPs  has  a  specific  category  assigned  to 


assess  the  actions  taken  subsequent  to  the  abandonment  of  the  MCR.  For  these  actions  a 
global  screening  value  of  1.0  is  assigned.  It  is  acknowledged  that  this  is  a  conservative 
estimate  and  that  more  detailed  analysis  should  be  performed  for  MCR  abandonment 
scenarios. 

The  scoping  method  introduced  in  NUREG-1921  offers  a  method  of  performing  a 
less  conservative,  but  still  simplified,  analysis  by  taking  into  account  how  specific  aspects 
of  the  fire  scenario  determine  operator  performance.  Using  the  decision-tree  format 
presented  in  NUREG-1921,  it  becomes  immediately  apparent  that  not  having  PORV 
indication  available  at  the  RSP  is  problematic.  For  our  scenario,  the  scoping  HRA 
decision  tree  simplifies  to  Figure  9-1. 


Necessary  Cues  for  Required  Actions  Protected 


Necessary  Cues  for  Required  Actions  NOT  Protected 


<  HEP=1.0 


Figure  9-1:  Scoping  HRA  Analysis  for  MCR  Abandonment  Scenario 


Individual  plants  that  show  susceptibility  to  MCR  abandonment  scenarios  may  be 
able  to  obtain  less  conservative  results  using  one  of  the  detailed  analysis  techniques 
presented  in  NUREG-1921.  However,  for  this  scenario,  the  absence  of  a  cue  to  alert  the 
operators  that  a  PORV  is  open  and  the  further  complication  caused  by  not  having  the 
PORV  block  valve  controls  on  the  RSP  electrically  independent  of  the  MCR  benchboard 
make  it  difficult  to  justify  taking  credit  for  successful  operator  action  after  abandoning 


the  MCR. 
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The  results  presented  in  NUREG-1 150  do  not  offer  a  methodology  for  obtaining 
the  probability  that  operators  will  successfully  recover  the  plant  from  the  RSP,  which 
requires  closing  the  PORV  block  valve.  The  value  is  given  in  the  NUREG-1 150  analysis 
as  a  Maximum  Entropy  distribution  characterized  by  a  lower  bound  of  a  =  0.0074,  an 
upper  bound  of  b  =  0.74,  and  a  mean  of  p  =  .074  such  that  the  probability  density 
function  (pdf)  is  given  by  [21]; 

f{e\a,b,M)  =  J, 

Where  p{/3  ^  0)  satisfies: 

1 

resulting  in  the  distribution  given  in  Eigure  9-2. 


Probability  of  Successful  Operator  Action  from  the 
RSP  Given  in  NUREG-  1150 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

Rop 


Figure  9-2:  Probability  of  Successful  Operator  Actiou  from  tbe  RSP  Giveu  iu  NUREG-1150 
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Although  not  taking  credit  for  operator  aetions  after  forced  abandonment  of  the 
MCR  revises  the  total  CDF  estimate  for  this  scenario  upward,  the  eombined  effeets  of 
reduced  fire  ignition  frequency  and  probability  of  MCR  abandonment  cause  our  updated 
distribution  of  CDF  for  this  scenario  to  be  signifieantly  less  than  the  values  predicted  in 
the  NUREG-1 150  analysis.  These  results,  along  with  a  discussion  of  the  most  important 
input  parameters,  will  be  presented  in  Seetion  X. 
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X,  CONCLUSIONS  AND  SENSITIVITY  ANALYSIS 


The  final  steps  in  our  analysis  are  to  combine  the  uncertainties  from  each  step  in 
our  scenario  to  develop  an  updated  CDF  distribution  that  takes  into  account  the 
interdependencies  between  the  input  parameters  and  to  identify  the  factors  that  have  the 
strongest  influence  on  the  final  results. 

In  order  to  obtain  an  overall  CDF  distribution  it  is  necessary  to  sample  from  the 
distributions  of  the  individual  factors  obtained  in  the  previous  sections.  From  Section  II, 
the  resulting  combined  CDF  is  calculated  as: 


Where: 

CDF  =  The  fire-induced  core  damage  frequency  for  the  MCR. 

=  The  frequency  of  MCR  fires. 

=  The  area  ratio  of  benchboard  1-1  to  total  cabinet  area  within  the  MCR. 
=  The  probability  that  operators  will  not  successfully  extinguish  the  fire 
before  forced  abandonment  of  the  MCR. 

=  The  probability  that  operators  will  unsuccessfully  recover  the  plant 

from  the  RSP.  To  successfully  recover  the  plant,  the  operator  must 
shut  the  PORV  block  valve,  despite  not  having  indication  on  the  RSP 
that  the  PORV  has  lifted. 


The  resulting  combined  CDF  distribution  is  given  in  Figure  10-1. 
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Figure  10-1:  Comparison  of  CDF  Distributions 


Despite  the  inability  to  take  eredit  for  operator  aetions  after  foreed  abandonment 
of  the  MCR,  each  of  the  current  methods  used  in  this  study  shows  that  a  reduction  by  a 
factor  of  approximately  two  in  the  mean  value  of  the  total  CDF  calculated  in  NUREG- 
1 150  is  expected.  This  is  shown  in  Table  10-1. 


CORE  DAMAGE  FREQUENCY  (1/yr) 

METHOD 

MEAN 

95th  PERCENTILE 

NUREG-1150 

1.58  X  10  ® 

1.98  X  10  ® 

PFS 

9.61  X  10'^ 

2.67  X  10  ® 

FMAG 

8.33  X  10'^ 

2.33  X  10  ® 

UMD 

6.88  X  10'^ 

1.92  X  10  ® 

Table  10-1:  Comparison  of  CDF  by  Method 
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Although  current  analysis  methods  indicate  an  overall  reduetion  in  the  mean  value 


of  eore  damage  frequeney,  the  differences  in  the  HEP  quantifieation  results  between 
NUREG-1 150  and  NUREG-1921  obscure  a  dramatic  decrease  in  the  expeeted  probability 
of  a  fire  in  benehboard  1-1  foreing  MCR  abandonment.  By  negleeting  the  Rop  term  our 
eore  damage  frequeney  equation  reduees  to  the  probability  of  operator  abandonment  of 
the  MCR  due  to  a  fire  in  benehboard  1-1,  sueh  that: 

^^{abandonment)  =  ^mcr  ^  fa  ^  fr 

Where: 

=  The  frequeney  of  MCR  fires. 

=  The  area  ratio  of  benehboard  1-1  to  total  eabinet  area  within  the 
MCR. 

=  The  probability  that  operators  will  not  successfully  extinguish  the  fire 
before  foreed  abandonment  of  the  MCR. 

Eaeh  of  the  eurrent  methods  used  to  ealeulate  the  probability  of  foreed  MCR 
abandonment  from  a  fire  in  benehboard  1-1  results  in  a  reduction  by  a  factor  of 
approximately  20  in  the  mean  value  of  the  probability  of  forced  MCR  abandonment 
calculated  in  NUREG-1 150.  This  is  shown  in  Table  10-2  and  Eigure  10-2. 


MCR  ABANDONMENT  PROBABILITY 

METHOD 

MEAN 

95th  PERCENTILE 

NUREG-1 150 

1.76  X  10  ® 

7.36  X  10  ® 

PFS 

9.61  X  10'^ 

2.67  X  10  ® 

FMAG 

8.33  X  10'^ 

2.33  X  10  ® 

UMD 

6.88  X  10'^ 

1.92  X  10  ® 

Table  10-2:  Comparison  of  MCR  Abandonment  Probability  by  Method 
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PROBABILITY  OF  MCR  ABANDONMENT  FROM  A 
BENCHBOARD  FIRE 
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Figure  10-2:  Comparison  of  MCR  Abandonment  Distributions 
Although  the  frequency  of  fires  causing  MCR  abandonment  is  much  lower  than 

previously  assessed,  it  is  still  a  significant  contributor  to  total  CDF  and  would  not  be 

eliminated  in  the  screening  process  of  a  fire  PRA. 

In  order  to  determine  which  input  parameters  have  the  strongest  influence  on  the 

results  of  our  fire  model,  PFS  uses  the  Spearman  Rank-order  Correlation  Coefficient 

(RCC)  [15]  to  assess  the  sensitivity  of  an  output  value,  Y,  to  an  input  parameter  X.  The 

RCC  is  defined  as: 


RCC 


=  1  - 


df 

n(n^  -  1) 


Where: 

di  =  X,  -  y,  =  The  difference  between  ranks  of  each  observation. 
n  =  The  number  of  data  pairs. 
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RCC  measures  the  degree  of  monotonieity  between  input  parameters  and  the 
chosen  output,  increasing  in  magnitude  to  a  value  of  unity  as  input  and  output  values 
approach  perfect  monotone  functions  of  each  other.  A  positive  value  of  RCC  indicates 
that  as  X  increases,  Y  tends  to  increase.  A  negative  value  of  RCC  indicates  that  as  X 
increases,  Y  tends  to  decrease. 

RCC  replaces  raw  scores  with  their  associated  ranks  to  reduce  the  effects  of 
nonlinear  data  [27].  Therefore,  RCC  is  independent  of  the  distribution  of  the  input 
parameters  and  allows  the  simultaneous  identification  of  both  modeling  parameters  and 
MCR  properties  that  have  the  strongest  influence  on  the  peak  HGL  temperature  achieved 
during  our  scenario  [10]. 

The  sensitivity  of  peak  HGL  temperature  to  the  input  parameters  varied  in  Table 
6-1  is  shown  in  Figure  10-3.  As  expected  [9],  HRR  has  the  most  direct  effect  on  peak 
HGL  temperature,  with  HRR  growth  time  and  ventilation  rate  also  having  significant 
effects.  Thermal  properties  of  the  MCR  are  minor  factors. 


Rank  Order  Correlation  for  Peak  HGL  Temperature 
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Figure  10-3:  Rank  Order  Correlation  Coefficient  for  Peak  HGL  Temperature 
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HRR  and  HRR  growth  time  are  modeling  parameters  that  are  out  of  the  operator’s 


control,  whereas  ventilation  can  be  easily  varied.  Figure  10-4  illustrates  how  varying 
MCR  ventilation  rate  for  the  FMSNL  21  test  case  modeled  in  Section  IV  affects  the  FIGL 


temperature. 


HGL  Temperature  for  Different  Ventilation  Rates 
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Figure  10-4:  FMSNL  21  HGL  Temperatures  as  a  Fuuctiou  of  Veutilatiou  Rate 
From  Figure  10-4  it  is  apparent  that  a  MCR  abandonment  condition  due  to 
exceeding  the  peak  HGL  temperature  threshold  will  not  be  predicted  by  CFAST  if  the 
ventilation  rate  is  increased  from  one  to  three  room  changes  per  hour  in  the  FMSNL  21 
test  case. 

The  FMSNL  21  test  case  uses  a  HRR  of  470  kW  that  is  greater  than  92%  of 
expected  fires.  To  find  a  ventilation  rate  that  would  more  nearly  preclude  abandonment 
due  to  peak  HGL  temperature  in  a  MCR  identical  to  the  FMSNL  21  test  case,  a  limiting 
case  simulation  was  conducted  with  a  HRR  of  702  kW  (98*'^  percentile)  and  a  minimal 
HRR  growth  time.  In  this  extreme  case,  a  ventilation  rate  of  10  room  changes  per  hour 


69 


was  sufficient  to  prevent  the  predicted  peak  HGL  temperature  from  reaching  the 
abandonment  threshold  of  95  °C. 


In  addition  to  our  tire  model,  our  analysis  is  sensitive  to  other  factors.  From  the 
CDF  equation  introduced  in  Section  II  it  is  apparent  that  refinements  of  MCR  fire 
ignition  frequency  estimates  would  have  a  linear  effect  on  the  estimate  CDF.  Also,  the 
effect  of  increasing  the  rate  at  which  fires  are  extinguished  in  the  MCR  (^mcr)  would 
cause  the  mean  probability  of  suppressing  a  fire  that  would  otherwise  cause  MCR 
abandonment  to  increase.  As  an  example,  the  effect  of  doubling  the  rate  at  which  fires 
are  extinguished  increases  the  mean  suppression  probability  by  a  factor  of  1 .29.  This  is 
shown  in  Figure  10-5. 
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Figure  10-5:  Effect  of  Doubling  Fire  Suppression  Rate 

Refinements  in  fire  ignition  frequency  and  fire  suppression  estimates  are  expected 

to  revise  the  severity  of  this  casualty  downward  as  measures  to  prevent  and  suppress  fires 
improve  with  time  due  to  increased  training  and  knowledge  sharing  practices. 
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The  scoping  method  of  quantifying  the  HEP  introduced  in  NUREG-1921[12]  is 
an  adequate  and  conservative  tool  for  our  analysis.  For  our  scenario,  the  design  of  the 
RSP  is  the  primary  factor  leading  to  our  conservative  HEP  estimate  of  1.0.  If  proper 
indication  were  provided  to  operators  once  they  were  forced  to  abandon  the  MCR,  and  if 
their  required  actions  did  not  involve  operation  of  equipment  not  electrically  independent 
of  the  fire,  it  is  likely  that  the  HEP  estimate  would  be  reduced  by  a  factor  of  2.  However, 
this  analysis  would  involve  a  review  of  plant  specific  procedures  and  time  margins  to 
core  damage. 

The  results  of  our  case  study  indicate  that  for  existing  plants  the  one  controllable 
factor  available  to  an  operator  to  mitigate  the  probability  of  a  MCR  abandonment 
scenario  is  the  MCR  ventilation  rate. 

For  plants  yet  to  be  constructed,  the  effects  of  a  MCR  casualty  could  be 
substantially  reduced  through  an  improved  design  of  the  RSP  that  provides  the  operators 
with  the  necessary  cues  and  independent  circuitry  to  take  actions  to  terminate  this 
casualty  before  core  damage  occurs. 

The  results  of  this  study  are  heavily  dependent  on  the  distribution  of  HRR.  The 
research  into  HRR  distributions  that  is  currently  being  conducted  will  help  to  provide 
regulators  with  the  necessary  tools  to  fully  assess  the  contribution  to  core  damage 
frequency  from  fires  initiated  from  within  the  MCR. 

It  is  important  to  note  that  this  work  is  a  limited  study  meant  to  apply  to  a  specific 
scenario  and  may  not  be  applicable  to  all  scenarios,  such  as  a  complete  loss  of  power  or  a 
fire  in  a  vertical  cabinet.  PRA  practitioners  are  encouraged  to  take  note  of  the  narrow 
range  of  applicability  of  the  methods  and  results  presented  here. 
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APPENDIX  A;  FMSNL  INPUT  FILE 


VERSN, 6, FM  Test  21 

1 1 

!! Environmental  Keywords 

1 1 

TIMES,  1800, -50,0,  10,  1 
EAMB, 288 . 15,101300,0 
TAMB, 288.15,101300,0,50 
C JET, WALLS 
CHEMI,10,393.15 

WIND, 0,10, 0.16 
1 1 

! ,  ! ,  Compartment  keywords 

COMPA, Compartment  1,18.3,12.2,6.1,0,0, 0, MariniteFM, Concrete FM, 
MariniteFM 
1 1 

!  !vent  keywords 

1 1 

WENT  ,2,1,1.08,2,1 

MVENT, 2, 1, 1,H, 4. 9,0. 66, H, 4. 9, 0.66, 0.38, 200, 300,1 
1 1 

!!fire  keywords 

I , 

OBJECT, FMSNL_21, 1, 12, 6. 1,0, 1,1, 0,0, 0,1 
1 1 

,!,! target  and  detector  keywords 

DETECT, 1,1, 347. 04, 3. 05, 6. 1,5. 98,100,0, 7E-05 

DETECT, 1,1,347.04,15.25,6.1,5. 98,100,0, 7E-05 

TARGET, 1,12,6.1,5.66,0,0,-1,TC, IMPLICIT, ODE 

1 1 

! ! misc .  stuff 

1 1 

THRMF, thermalfmsnl 
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APPENDIX  B:  FMSNL  21  INPUT  TO  CFAST-SCREEN  VIEW 


1 .  Simulation  Environment 


File  Run!  Tools  View  Help 

Simulation  Environment  j  Comoartment  Geometrv  1  Horizontal  Flow  Vents  |  Vertical  Flow  Vents  |  Mechanical  Flow  Vents)  Fires  |  Detection /SuDoression )  Taraets )  Surface  Connections ) 
-  “Ambient  Conditions - 


Title  |FM  T6st21 


“Simulation  Times 

Simulation  Time: 

)TdOOs 

T©cl  Output  interval: 

|50s 

Binary  Output  Interval 

)0s 

Spreadsheet  Output  Interval: 

jl  s 

Smokeview  Output  interval 

)10S 

rThermel  Properties  File  “ 
jthermalfmsnl.csv 


■J 


“Intenor 

Temperature 

|15‘C 

Elevation 

)0m 

Pressure 

(101300  Pa 

Humidity. 

)50% 

Temperature.  |l5  'C 


Elevation:  [Orr 


Pressure  |l01300Pa 

Wind  Speed  jO  m/s 
Scale  Height  jlO  m 


No  Errors  orWamings 


-  Input  File  Syntax  Check  0 


2.  Compartment  Geometry 


File  Run!  Tools  View  Help 

Simulation  Environment  Compartment  Geometry  ]  Horizontal  Flow  Vents  |  Vertical  Flow  Vents  |  Mechanical  Flow  Vents)  Fires  |  Detecton/ Suppression )  Taraets )  Surface  Connections ) 


Duplicate 


Move  Up 


“Compartment  1  (of  1)  “ 

“Geometry - 


Compartment  Name:  |Compartment  t 


“Flow  Characteristics 


Width  (X): 

Depth  (Y): 

Height  (Z): 

ll8.3m 

Position,  X: 

|0  m 

|l2.2m 

Y- 

|0  m 

|6.1  m 

Z 

|0  m 

Ceiling: 

(Marinite  FMSNL 

“3 

Walls 

(Marinite  FMSNL 

“3 

Floor 

)Concrete  FMSNL 

“3 
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3.  Horizontal  Flow  Vents 


File  Run!  Tools  View  Help 

Simulation  Environment  I  ComoartmentGeome&v  Horizontal  Flow  Vents  j  Vertical  Flow  Vents]  Mechanical  Flow  Vents)  Fires  )  Detecton/SuoDression)  Taroets)  Surface  Connections) 


Num  I  First  Compartment  |  Offset  1  [Second  Compartmgij  Offset  2  [  Sill  |  Soffit  |  Width  |  Wind  Initial  Opeij~ 


Sill  r 

Soffit  r 


Width:  r 


Duplicate  )  Move  Up 


First  Compartment 

“Second  Compartment 

1  d 

1  d 

Vent  Offs^:  |0  m 

Vent  Offset:  )0m 

Initial  Opening  Fraction  |l~ 
Change  Fraction  At:  | 
Final  Opening  Fraction'  | 


Wind  Angle 

Face  [Front” 


“3 


4.  Vertical  Flow  Vents 


File  Run!  Tools  View  Help 

Simulaton  Environment)  Comoartment  Geometry  j  Horizontal  Flow  Vents  Vertical  Flow  Vents  )  Mechanical  Flow  Vents)  Fires  )  Detection /Suooression )  Taroets )  Surface  Connections ) 


Add  I 


Duplicate  | 


“Vent  1  (of  1)  Geometry 


Top  Compartment 

“Bottom  Compartment 

)Outside 

)  Compartment  1 

Initial  Opening  Fraction:  |i 

Cross-Sectional  Area:  )l  .08m^ 

Change  Fraction  At:  (oT 

Shape  [square  ^ 

Final  Opening  Fraction:  [i 


Save  [  Run  [  View 

No  Errors 
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5.  Mechanical  Flow  Vents 


File  Run!  Tools  View  Help 

Simulator!  Environment  I  Comoartment  Geome&v  1  Horizontal  Flow  Vents  |  Vertcal  Flow  Vents  Mechanical  Flow  Vents  [Fires  j  Delecton/SuDoressionl  Taraetsj  Surface  Connections] 


“Vent  1  (of  1)  Geometry 


Flow  Rate 

|0.38  m'S/s 

Initial  Opening  Fraction:  P 

Begin  Dropoff  Af 

[200  Pa 

Change  Fraction  At:  jos 

Zero  Flow  At: 

[300  Pa 

Final  Opening  Fraction:  h 

Save  I  Run  j  View 


6.  Fires 


File  Run!  Tools  View  Help 


Simulation  Environment!  Comoartment  Geomebv  |  Horizontal  Flow  Vents  |  Vertcal  Flow  Vents  |  Mechanical  Flow  Vents  Fires  ]  Detection /SuDoression  j  Taraets  j  Surface  Connections  j 


Ceiling  Jrt:  jCeiling  &  Walls" 
Lower  Oxygen  Limit:  |l0% 


Gaseous  Ignition 

Temperature  '  ^ 


Duplicate 


Compartment.  jCompartmentl 


Type 

jConstrained  ^ 

Position.  X  |l2m 

Position  Y:  |6.1  m 

Position  Z  |0  m 

Ignition  Critwion  jTime 

Normal,  X; 

l» 

Normal.  Y:  |o 

Normal.  Z'  |l 

Plume.  jMcCaffrey 

Ignition  Value:  jo  s 

“Fire  Object - 

Fire  Object:  |FMSNL_21  ^  _ Edit  j 

Material:  Methane,  a  transparent  gas  (CH4) 

Length  1  m 

Width:  1  m 

Thickness:  0.25  m 

MolarMass  0.1002kg/mol 

Total  Mass  10000  kg 

Heal  of  Combustion:  45000  kJ/kg 

H^t  of  Gasification:  0  kJ/kg 

Vcrfitilization  Temperature:  22  *C 

Radiative  Fraction:  0.35 


FMSNL_21  HRR 


No  Btws 
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7.  Detection/Suppression 


Simulation  Environment  I  ComoartmentGeome&vl  Horizontal  Flow  Vents  |  Vertical  Flow  Vents  |  Mechanical  Flow  Vents)  Fires  Detection /Suppression  llaraets)  Surface  Connections) 


)  Num  Compartment 

Type 

1  X  Position 

Y  Position  . 

Z  Position 

Activation 

RTI  1 

Spray  Density 

_ 1  * 

1  I  CompartrT>ent  1 

Smoke 

1  305 

61 

598 

73  89001 

100 

7E-05 

H 

2  Compartment  1 

Smoke 

15.25 

61 

5.98 

73  89001 

100 

7E-05 

Add  I  Duplicate  |  Move  Up  |  Move  Down  | 


-Alarm  1  (af2)  - 


Type.  )Smoke  Alarm  ^  Compartment:  )  Compartment  1 


Width  (X): 

j3.05m 

Depth  (Y); 

j6.1  m 

Height  (Z): 

|5.98m 

■r)  Activation  Temperature  )73.89001  ‘C 


RTI:  )l00(ms)'^5 


Spray  Det^ity  )7E-05m/s 


No  Biws 


8.  Targets 

File  Run!  Tools  View  Help 

Simulaton  Environment)  Compartment  Geometry  |  Horizontal  Flow  Vents  |  Vertical  Flow  Vents  |  Mechanical  Flow  Vents)  Fires  )  Detection /Suppression  Targets  j  Surface  Connections ) 


Add  I  Duplicate  |  Move  Up  |  Move  Down  |  Remove  | 


"Target  1  (of  1)  Geometry - 

Compartment  )Comparmient1  ^ 


"Position 

"Normal  Vector  Points  To 

Matenal:  )Thermocouple  (small  steel  target  far 

Width  (X):  |l2m 

Depth  (Y):  |6.1  m 

Height  (Z):  |5.66  m 

)User  Specified  ^ 

Wxtth  (X):  )0 

Depth  (Y);  )0 

Method;  jimplicit 

Type:  jThermally  Thin  J 

Height  (Z):  )-1 

Save  I  Run  |  View 
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9.  Surface  Connections 


Fite  Run!  Tools  View  Help 

Simulation  Environment  I  ComoartmentGeome&vl  Horizontal  Flow  Vents  |  Vertical  Flow  Vents  ]  Mechanical  Flow  Vents)  Fires  )  Detection /SuDoression)  Taraets  Surfece  Connections  | 


■Vertical  Connections  “ 


HeotTransferConnectton  1  ' 

Top  Compartment  | 

Bottom  Compartment  | 


U 
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APPENDIX  C:  FMSNL  THERMAL  PROPERTIES  INPUT  FILE 


Short  Name 

Conductivity 

Specific 

Heat 

Density 

Thickness 

Emissivity 

HCi 

Coefficients 

Long  Name 

METHANE 

0.07 

1090 

930 

0.0127 

0.04 

0 

Methane,  a  transparent  gas  (CH4) 

MineralBE2 

0.2 

150 

500 

0.05 

0.95 

0 

Minerai  Wooi  BE2 

SteelBE2 

54 

425 

7850 

0.001 

0.95 

0 

Steei  iCFMP  BE2 

ConcreteBE2 

2 

900 

2300 

0.15 

0.95 

0 

Concrete  iCFMP  BE2 

MARIBE3 

0.12 

1250 

737 

0.0254 

0.8 

0 

Marinite  iCFMP  BE3  (2  1/2  in  iayers) 

GYPBE3 

0.16 

900 

790 

0.0254 

0.9 

0 

Gypsum  iCFMP  BE3  (2  1/2  in  iayers) 

XLP  C  BE3 

0.21 

1560 

1375 

0.01 

0.95 

0 

FXPE  Cabie  iCFMP  BE3 

PVC  C  BE3 

0.147 

1469 

1380 

0.01 

0.95 

0 

PVC  Siab  iCFMP  BE3 

XLP  P  BE3 

0.21 

1560 

1375 

0.191 

0.95 

0 

F  XPE  Cabie  iCFMP  BE3 

PVC_P_BE3 

0.147 

1469 

1380 

0.191 

0.95 

0 

PVC  Siab  iCFMP  BE3 

SteelBE4 

44.5 

480 

7743 

0.02 

0.95 

0 

Steei  iCFMP  BE4 

ConcreteBE4 

2.1 

880 

2400 

0.25 

0.95 

0 

Concrete  iCFMP  BE4  and  BE5 

LiteConcBE4 

0.11 

1350 

420 

0.3 

0.95 

0 

Lightweight  Concrete  ICFMP  BE4/BE5 

PVC  P  BE4 

0.134 

1586 

1380 

0.015 

0.8 

0 

PVC  Power  Cabie  BE4 

PVC  C  BE4 

0.134 

1586 

1380 

0.007 

0.8 

0 

PVC  Controi  Cabie  BE4 

MariniteFM 

0.12 

1250 

720 

0.025 

0.95 

0 

Marinite  FMSNL 

ConcreteFM 

1.8 

1040 

2280 

0.15 

0.95 

0 

Concrete  FMSNL 

FireBrickNBS 

0.36 

1040 

750 

0.113 

0.8 

0 

Fire  Brick  NBS 

CeramicNBS 

0.09 

1040 

128 

0.05 

0.97 

0 

Ceramic  Fiber  NBS 

MariniteNBS 

0.12 

1250 

720 

0.0127 

0.83 

0 

Marinite  NBS 

GypsumNBS 

0.17 

1090 

930 

0.0127 

0.95 

0 

Gypsum  NBS 

ConcreteNBS 

1.8 

1040 

2280 

0.102 

0.95 

0 

Concrete  NBS 

TC 

54 

425 

7850 

0.001 

0.95 

0 

Thermocoupie  (smaii  steei  target  for 
piume  temp) 
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APPENDIX  D:  WINBUGS  INPUT  FILE  FOR  HGL  DATA  (BE3  ONLY) 


bm~dunif(-10,10) 
sm~  dunif(0,10) 
taum<- l/pow(sm,2) 
pe<-0.08 

be<-(log(  1  +pe)+log(  1  -pe))/2 
se<-(log(l+pe)-log(l-pe))/(2*  1.95996398454005) 

bt<-bm-be 

st<-sqrt(pow(sm,2)+pow(se,2)) 

C  <-  1000 

for(  i  in  1  ;  N  ) 

{ zeros  [i]  <-  0 

L[i]  <-  pow(exp(-0.5*pow((log(x[i,2]/x[i,l])- 
bt)/st,2))/(sqrt(2*3.141592654)*st)/(x[i,2]/x[i,l]),x[i,3]) 
ghr[i]  <-(-l)Mog(L[i])  +  C 
zeros[i]  ~  dpois(ghr[i])} 

fm~dlnorm(bm,taum) 

logfm<-log(fm) 

for(  j  in  1  ;  9  ) 

{yyD']<-xxDTfm 

P2 1 8  [)]<- 1  -phi((log(2 1 8)-bm-log(xx[j]))/sm) 
P330[j]<-l-phi((log(330)-bm-log(xx[j])ysm)} 
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APPENDIX  E:  WINBUGS  INPUT  FILE  FOR  HGL  DATA  fBE3  +  FMSNL  21/22) 


bm~dunif(-10,10) 
sm~  dunif(0,10) 
taum<- l/pow(sm,2) 
pe<-0.08 

be<-(log(  1  +pe)+log(  1  -pe))/2 
se<-(log(l+pe)-log(l-pe))/(2*  1.95996398454005) 

bt<-bm-be 

st<-sqrt(pow(sm,2)+pow(se,2)) 

C  <-  1000 

for(  i  in  1  ;  N  ) 

{ zeros  [i]  <-  0 

L[i]  <-  pow(exp(-0.5*pow((log(x[i,2]/x[i,l])- 
bt)/st,2))/(sqrt(2*3.141592654)*st)/(x[i,2]/x[i,l]),x[i,3]) 
ghr[i]  <-(-l)Mog(L[i])  +  C 
zeros[i]  ~  dpois(ghr[i])} 

fm~dlnorm(bm,taum) 

logfm<-log(fm) 

for(  j  in  1  ;  9  ) 

{yyD']<-xxDTfm 

P2 1 8  [)]<- 1  -phi((log(2 1 8)-bm-log(xx[j]))/sm) 
P330[j]<-l-phi((log(330)-bm-log(xx[j])ysm)} 
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APPENDIX  F:  FMSNL  21  INPUT  TO  PFS 


I.  CF AST  CONTROL 


CFAST  Model  and  Execution  Control 


PFS:  Probabilistic  Fire  Simulator  v4.0 
VTT  Technical  Research  Centre  of  Finland 


Fire 


Sprinkler/Detector  Control 


Type 

1 

0  =  constant,  1  =  t2  +  exp.  dr 

Detector  type 

1 

1=smoke,  2  =  heat 

2  =  Flame  front  velocity 

Sprinkler 

0 

0  =  detector,  1  =  sprinkle 

4  =  Cable  fire  (Mangs) 

RTI 

100 

(ms)''(1/2) 

RHR  Growth  time  (s) 

320 

101-114  Household  appl. 

Tact 

73.89001 

C 

RHR  limiter  (kW) 

470 

201  -205  Detector  test  fires 

Decay  rate  (s) 

10 

301-306  NPP  cabinets 

Decay  start  (s) 

1140 

MaxRHR  (kW) 

470.0 

CFAST  control 

Scaling  factor 

1 

Cfast  directory 

C:\NIST\cfast51 1 

Iteration  control 

xval 

nsteps 

Source  height  (m) 

0 

DTCHECK 

1.0E-09 

1001 

(default 

1.0E-09 

100) 

RHR/fuel  area 

470 

kW/m2 

Nrows 

101 

18|dt  output 

CallMode 

1 

0  =  normal  2  =  save 

Ambient 

1  =  debug 

Create  thermal  database 

Ambient  temp 

15 

C 

ThDbCreate 

TRUE 

(default  FALSE) 

Ambient  pres 

1.013 

bar 

If  thermal  data  is  randomly  sampled, 

choose  TRUE.  Otherwise  use  FALSE 

and  existing  thermal  data  base  is  used  (faster). 
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II.  CFAST  INPUT 


CFAST  Model  Input  Sheet 

Create  CFAST  Input  Data  File 

Current  Iteration 

0 

Running  directory 

C:\NIST\PFS 

1 

NROOMS 

DAT  file 

CRAB  O.dat 

HI  file 

CRAB 0.hi 

TARGET  ROOMS 

Thermal  DB 

Tunnel 

(Do  not  include  file  extension!) 

Thermal  DB  to  use 

CRAB  0 

1  2 

3 

First  target  room  contains  the  heating  target  (cable). 

Intervals 

Time  parameters 

(S) 

Total  time 

Print 

History 

Display 

Restart 

Maximum  RunTime 

TIMES 

1800 

180 

1 

0 

180|s 

Ambient  conditions 

Temp  (C) 

Pressure  (P  Elev.  (m) 

TAMB 

15 

101300 

0 

EAMB 

15 

101300 

0 

Room  # 

Room  geometry 

1 

SUM 

Hl/F 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Y 

WIDTH 

12.2 

X 

DEPTH 

18.3 

18.3 

Z 

HEIGH 

6.1 

Area 

223.26 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

223.26 

Volume 

1361.69 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1361.89 

Index  to  the  thermal  data  base. 

Ceiling  and  wall 

CEiLI 

3 

3 

3 

4 

4 

4 

4 

4 

4 

4 

materials 

WALLS 

3 

3 

3 

4 

4 

4 

4 

4 

4 

4 

Vent  connect 


HVENT 

HVW 


HVO  Soffit 

(m) 


Second 

Vent 


(m) 


Second 

Vent 


Window  breaking  temperature  Q 


TARGETS 
NTargets  Q 


1  These  default  targets  are  defined  to  have  the  same  position,  but  different  oriantation,  to  take  the  direction  dependence  into  account. 


Room#  X(m)  Y  (m)  Z(m)  Normal  X  Normal  Y  Normal  Z  Material#  Method  Eq.  Type 

1 

1  12.00  6.10  5.66  0  0  -1  10 

IMPLICIT  PDE  -  -• 

IMPLICIT  PDE 
IMPLICIT  PDE 
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DETECTORS 


1  2 

(1  for  smoke,  2  for  heat) 

DETECT 

Type  Room#  Tactiv(C)  X{m)  Y(m)  Z(m)  RTI  {m.s)1  Sprinkler  Spray  dens 

1 

1 

1 

73.89001 

3.05 

6.10 

5.98 

100 

0 

0 

2 

1 

1 

73.89001 

15.25 

6.10 

5.98 

100 

0 

0 

CAUTION!  Don't  use! 
spray  density  >  0.0046  gives  extiction, 
0.0046  nothing!! 


FANS  It  is  possible  to  define  fans  between  the  rooms,  but  not  ductwork 


Nfans 


1  MVOPN  1  =  Horizontal  (in  ceiling),  2  =  vertical  MVFAN 


1st  room  Orientatiion  Height  (m) 

2nd  room  Orientation  Height  (m) 

Area  (m2)  Pmin(Pa)  Pmax(Pa)  Flow(m3/s] 

1 

1 

1 

4.900 

2 

1 

4.900 

0.66 

0.00 

300 

0.38 

t  - r—y 

FIRE  SOURCE  PROPERTIES 


Chemistry 


CHEMI 


Molar  weig  Humidity 
(%) 

LOI 

{%) 

-DHc 

(J/kg) 

T(fuel,init)  T(ign,  gas)  Radiative 

K  K  fraction 

16 

50 

10  4.50E+07 

288.15  393.15  0.35 

Compartment  of  fire  origin 

LFBO  I  1| 


Fire  Type  _ 

LFBT  I  ^  1  =  unconstrained  fire,  2  =  constrained 


Fire  Position 

FPOS 


X(m)  Y{m)  Z(m) 
12|  5.1|  0 


Maximum  Fire  Area 

I  ]]m2 


Fire  height 

I  °i'" 


Ceiling  jet 

CJET 


OFF.  CEILING,  WALL,  ALL 

I 


Time  series _ Speciel  Specie2 


FTIME 

FMASS  FHIGH 

FAREA 

FQDOT 

OD  CO 

0 

0 

0 

0 

1  0 

0.01 

0 

1 

90 

0.001758 

0 

1  79101.56 

0.01 

0 

2 

180 

0.007031 

0 

1  316406.3 

0.01 

0 

3 

270 

0.010444 

0 

1  470000 

0.01 

0 

4 

360 

0.010444 

0 

1  470000 

0.01 

0 

5 

450 

0.010444 

0 

1  470000 

0.01 

0 

6 

540 

0.010444 

0 

1  470000 

0.01 

0 

7 

630 

0.010444 

0 

1  470000 

0.01 

0 

8 

720 

0.010444 

0 

1  470000 

0.01 

0 

9 

810 

0.010444 

0 

1  470000 

0.01 

0 

10 

900 

0.010444 

0 

1  470000 

0.01 

0 

11 

990 

0.010444 

0 

1  470000 

0.01 

0 

12 

1080 

0.010444 

0 

1  470000 

0.01 

0 

13 

1170 

0.001106 

0 

1  49787.07 

0.01 

0 

14 

1260 

1.37E-07 

0 

1  6.144212 

0.01 

0 

15 

1350 

1.69E-11 

0 

1  0.000758 

0.01 

0 

16 

1440 

2.08E-15 

0 

1  9.36E-08 

0.01 

0 

17 

1530 

2.57E-19 

0 

1  1.15E-11 

0.01 

0 

18 

1620 

3.17E-23 

0 

1  1.43E-15 

0.01 

0 

19 

1710 

3.91E-27 

0 

1  1.76E-19 

0.01 

0 

20 

1800 

4.82E-31 

0 

1  2.17E-23 

0.01 

0 

Possible  species 

HCN 

HCL 

HCR 

CT 

02 

OD 

CO 
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APPENDIX  G:  STUDENT  BIOGRAPHY 


Lieutenant  Mark  Minton  is  originally  from  Northern  California  and  enlisted  in  the 
United  States  Navy  after  graduating  from  Clayton  Valley  High  Sehool  in  1995.  After 
eompleting  the  Naval  Nuelear  Power  Sehool  enlisted  eurrieulum,  he  was  seleeted  for  the 
Nuelear  Enlisted  Commissioning  Program  and  ordered  to  Oregon  State  University,  where 
he  reeeived  a  Baehelor  of  Seienee  in  Nuelear  Engineering  in  2001 . 

After  being  eommissioned,  he  reported  to  USS  SHOUP  (DDG  86),  where  he 
served  as  Communieations  Offieer  and  aehieved  qualifieation  as  a  Surfaee  Warfare 
Offleer.  He  was  subsequently  ordered  to  return  to  Naval  Nuelear  Power  Sehool,  this  time 
to  eomplete  the  offieer  eurrieulum. 

Upon  eompletion  of  Nuelear  Power  training,  Eieutenant  Minton  reported  to  USS 
RONAED  REAGAN  (CVN  76)  where  he  served  in  Reaetor  Department  and  aehieved 
qualifieation  as  Nuelear  Engineering  Offieer.  He  was  subsequently  seleeted  for  lateral 
transfer  to  the  Engineering  Duty  Offieer  (Nuelear)  eommunity. 

In  2008,  Eieutenant  Minton  reeeived  a  Master  of  Engineering  Management  from 
Old  Dominion  University. 

Eieutenant  Minton  is  eurrently  pursuing  a  Nuelear  Engineer’s  Degree  and  a 
Master  of  Seienee  in  Nuelear  Seienee  in  Engineering  at  the  Massaehusetts  Institute  of 
Teehnology. 

Upon  eompletion  of  his  studies  at  the  Massaehusetts  Institute  of  Teehnology, 
Eieutenant  Minton  has  been  ordered  to  report  to  the  Supervisor  of  Shipbuilding,  Newport 
News,  Virginia,  where  he  will  be  assigned  duties  pertaining  to  the  design,  eonstruetion, 
eonversion,  and  overhaul  of  nuelear  powered  aireraft  earners. 
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